Polyakov loop in chiral quark models at finite temperature 
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We describe how the inclusion of the gluonic Polyakov loop incorporates large gauge invariance and 
drastically modifies finite temperature calculations in chiral quark models after color neutral states 
are singled out. This generates an effective theory of quarks and Polyakov loops as basic degrees of 
freedom. We find a strong suppression of finite temperature effects in hadronic observables triggered 
by approximate triality conservation (Polyakov cooling) , so that while the center symmetry breaking 
is exponentially small with the constituent quark mass, chiral symmetry restoration is exponentially 
small with the pion mass. To illustrate the point we compute some low energy observables at 
finite temperature and show that the finite temperature corrections to the low energy coefficients 
are Nc suppressed due to color average of the Polyakov loop. Our analysis also shows how the 
phenomenology of chiral quark models at finite temperature can be made compatible with the 
expectations of chiral perturbation theory. The implications for the simultaneous center symmetry 
breaking-chiral symmetry restoration phase transition are also discussed. 

PACS numbers: ll.lO.Wx 11.15.-q ll.lO.Jj 12.38.Lg 

Keywords: finite temperature; heat kernel expansion; effective action; gauge invariance; chiral quark models; 
chiral Lagrangian; Polyakov Loop 



I. INTRODUCTION 

The general belief that QCD undergoes a phase transi- 
tion to a quark-gluon plasma phase at high temperature 
has triggered a lot of activity both on the theoretical 
as well as on the experimental side. The original argu- 
ment put forward by Casher ij suggesting that confine- 
ment implies dynamical chiral symmetry breaking and 
hence that the chiral and deconfinement phase transi- 
tions take place simultaneously at least at zero chemical 
potential, has been pursued and so far confirmed in theo- 
retical studies on the lattice (2| ■ This also agrees with the 
phenomenological determinations of the vacuum energy 
density in the bag model, with an energy density differ- 
ence between the Wigner and Goldstone realizations of 
chiral symmetry. It has also been shown that in the large 
Nc limit with the temperature T kept fixed, if a chiral 
phase transition takes place it should be first order Q . 

The coupling of QCD distinctive order parameters 
at finite temperature to hadronic properties has been 
the subject of much attention over the recent past 
0, IH IE 0) S 13 mainly in connection with theoreti- 
cal expectations on the formation of quark-gluon plasma 
and the onset of deconfinement. Indeed, even if such 
a state of matter is produced in existing (RHIC, SPS 
and future (LHC) facilities, the states which are 
detected are hadrons created in a hot environment. Thus, 
it makes sense to study the properties of hadrons in a 
medium which can undergo a confinement-dcconfinement 
phase transition. For heavy masses, quarks become static 
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sources and there is a general consensus that the order 
parameter can be taken to be the Polyakov loop or ther- 
mal Wilson line ^3 where the breaking of the center 
symmetry signals the onset of deconfinement. Dynami- 
cal light quarks, however, break explicitly the center sym- 
metry and no criterion for deconfinement has been estab- 
lished yet I'^ QCD, there has been an increasing 
interest in developing effective actions for the Polyakov 
loop as a confincmcnt-deconfinement order parameter be- 
cause of their relevance in describing the phase transition 
from above the critical temperature [la lla flTl Il8| . 

On the other hand, in a hot medium, one also ex- 
pects that the spontaneously broken chiral symmetry 
is restored at some critical temperature. For this chi- 
ral phase transition the quark condensate is commonly 
adopted as the relevant order parameter. The melting 
of the chiral quark condensate has been observed on the 
lattice 01 , is suggested by chiral perturbation theory ex- 
trapolations [T9I I20I I and is numerically reproduced in 
chiral quark models before pi, 22] and after inclusion of 
pion corrections "23] (for a review see e.g. Ref. |2^). 

Where theory has most problems is precisely in the 
interesting intermediate temperature regime around the 
phase transition, because both the lightest Goldstone 
particles and the Polyakov loop degrees of freedom should 
play a role, if they coexist. Up to now it is uncertain 
how the corresponding states couple to each other from 
a fundamental QCD viewpoint, hence some modeling is 
required. Based on previous works 0, H, and to com- 
ply with chiral symmetry it seems natural to couple chiral 
quark models and the Polyakov loop in a minimal way 
as an effective space dependent color chemical potential. 
The work in Ref. accounts for a crossover between 
the restoration of chiral symmetry and the spontaneous 
breaking of the center symmetry, reproducing qualita- 
tively the features observed on the lattice simulations 
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and which find a natural explanation in terms of dimen- 
sion two condensates . In this regard we want to argue 
below that the special role played by the gauge symmetry 
at finite temperature actually requires this coupling, and 
elaborate on the consequences of it when the quantum 
gluon effects are considered. 

The organization of the paper is as follows. We review 
some facts on large gauge symmetry at finite temperature 
in Sect, which are put into the context of chiral quark 
models. Next, we address the problem suffered by chiral 
quark models at finite temperature in Section IlIII where 
we argue that the origin of the difficulty is related to a 
defective treatment of the large gauge symmetry at finite 
temperature. Thus, to comply with gauge invariance at 
finite temperature one should at least couple the quarks 
to the Aq gluon field. We do this in Section Hvl This is 
equivalent to make the replacement 

do^do + iAo (1.1) 

which corresponds to an x-dependent chemical poten- 
tial coupling in the color fundamental representation. 
Obviously, this coupling introduces a color source into 
the problem for a fixed Aq field. In order to project 
onto the color neutral states we integrate over the Ao 
field, in a gauge invariant manner. In Section we de- 
scribe the consequences of such a coupling and projec- 
tion in chiral quark models for a variety of observables 
at the one quark loop approximation. Actually, as we 
will show, there is an accidental Z{Nc) symmetry in the 
model which generates a triality (super)selection rule at 
this level of approximation, from which a strong ther- 
mal suppression, 

0(^e-N.M/T-^ follows in the quenched 
approximation. This puts some doubts on whether chi- 
ral quark models do predict a chiral phase transition at 
realistic tem per atures as we advanced in previous com- 
munications [23, 123 ■ Corrections beyond one quark loop 
are discussed in Section lVll where we see that the suppres- 
sion at low temperatures actually becomes ©(e"™"/-^), 
very much along the expectations of Chiral Perturbation 
Theory (ChPT) Gluonic corrections and local cor- 

rections in the Polyakov loop are also analyzed in this 
section. In view of our discussions we illustrate in Sec- 
tion the situation with schematic dynamical calcula- 
tions involving quantum and local Polyakov loops in the 
unqucnched theory as compared to lattice studies. In 
Section rVIIII we extend these calculations to the region 
around the phase transition. Finally, in Section IIXI we 
summarize our points and draw our main conclusions. 

II. GAUGE INVARIANCE OF CHIRAL QUARK 
MODELS AT FINITE TEMPERATURE AND THE 
POLYAKOV LOOP 

In this section we review some relevant and naively 
very disparate concepts of gauge symmetry at finite tem- 
perature. Sect. Ill and the center symmetry in gluody- 
namics. Sect. IIIBI as well as the standard chiral quark 



models. Sect. Ill CI in order to fix our notation for the rest 
of the paper. Both subjects are well known on their own, 
although rarely discussed simultaneously, and the reader 
familiar with any of them may skip the corresponding 
subsections. Advancing the result of subsequent discus- 
sions made in latter sections, the basic Polyakov Chiral 
Quark Model is first introduced in Sect. Ill Dl The con- 
flict between both large gauge symmetry and chiral quark 
models is discussed in Sect. lIIII The solution to the prob- 
lem is elaborated in Sect. II VI where the coupling of the 
Polyakov loop to chiral quark models is motivated. 



A. Large gauge symmetry 

One of the most striking features of a gauge theory 
like QCD at finite temperatures is the non-perturbative 
manifestation of the non Abelian gauge symmetry. In- 
deed, in the Matsubara formalism of quantum field the- 
ory at finite temperature the space-time becomes a topo- 
logical cylinder: one introduces a compactified Euclidean 
imaginary time 30] and evaluates the path integral sub- 
jecting the fields to periodic or antiperiodic boundary 
conditions for bosons and fermions respectively in the 
imaginary time interval /3 = 1/T where T is the tem- 
perature. We use the Euclidean notation X4 — ixQ and 
A4{x,X4) = jAo(x, xq). Thus, only periodic gauge trans- 
formations, g(x,X4) — g{x,X4 + P), are acceptable since 
the quark and gluon fields are stable under these trans- 
formations. In the Polyakov gauge ^4^4 = with A4 a 
diagonal traceless Nc x Nc matrix, one has for the gauge 
SU{Nc) group element 

g{x4) = diag(e^2...„,T) (3.1) 

(X^j^i "-i ~ 0) the following gauge transformation on the 
A4 component of the gluon field 

A4^ A4 + 27rTdiag(nj) . (2.2) 

Thus, in this particular gauge, gauge invariance manifests 
as the periodicity in the A4 gluon field. This property 
is crucial and prevents from the very beginning from the 
use of a perturbative expansion in the gluon field, A4, 
at finite temperature. This large gauge symmetry^ can 
be properly accounted for by considering the Polyakov 
loop or untraced Wilson line as an independent degree of 
freedom, 

nXi + l/T 

n{x)=Texpi dx'4A4{x,x4) (2.3) 



Technically speaking the transformations 12. li may not be large 
in the topological sense (i.e., homotopically non trivial). This 
depends on the t opo logy of the spatial manifold as well as on 
the gauge group j2lj. They are topologically large within the 
Polyakov gauge. 
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where T indicates the Euclidean temporal ordering ope- 
rator and A4 the gluon field. Under a general periodic 
gauge transformation one gets 

n{x) ^ g{x)n{x)gHx) . (2.4) 

In the Polyakov gauge, which we assume from now on, fl 
becomes 

n{x) = e*-^*'^)/^ (2.5) 

and so it is invariant under the set of transformations 
(12.1(1 . The failure of perturbation theory at finite tem- 
perature in a gauge theory has generated lot of discus- 
sion in the past mainly in connection with topological 
aspects, Chern-Simons terms, anomalies, etc. In the 
case of the topological Chern-Simons term radiatively in- 
duced by fermions in 2 + 1 dimensions |32j it was puzzling 
to find, in the perturbative treatment, that the Chern- 
Simons quantization condition j33l | was violated at finite 
temperature [^EEl- It was subsequently shown that, 
within a non perturbative treatment, no contradiction 
arises )3Q| . In j37j .38J it was shown that a derivative ex- 
pansion approach, suitably defined at finite temperature, 
was appropriate to deal with this problem. We will use 
this approach in the present work. 

B. Center symmetry in gluodynamics 

In pure gluodynamics at finite temperature one can 
use the center of the gauge group to extend the periodic 
transformations to aperiodic ones |39l | , 

9(.S,^)^zgix,0), (2.6) 

so that z is an element of 1j{Nc)- An example of such a 
transformation (with z = e'^'^/^c^ jjj ^j^g Polyakov gauge 
is given by 

3(0:4) =diag(e*2""^"^^/^=), 

m = 1 - Nc, nj>2 = 1 (2.7) 

and the gauge transformation on the A4 component of 
the gluon field is 

2'kT 

Aa^ + — diag(nj) . (2.8) 

Under these transformations both gluonic action, mea- 
sure and boundary conditions are invariant. The 
Polyakov loop, however, transforms as the fundamental 
representation of the Z,{Nc) group, i.e. CI zfl, yielding 
(fi) = and hence = 0. More generally, in the 

center symmetric or confining phase 

(17") =0 for n^kNc, keZ. (2.9) 

Actually, this center symmetry is spontaneously bro- 
ken above a critical temperature, Td « 270 MeV for 



iVc = 3 The antiperiodic quark field boundary 

conditions are not preserved under non trivial center 
transformations since g(x, l/T) — > g{x, 1/T)q{x, l/T) = 
—zg{x, 0)g(x, 0) instead of —g{x, 0)q{x, 0). A direct con- 
sequence of such property is the vanishing of contribu- 
tions to the quark bilinear of the form^ 

{q{n/T)q{Q)) = for kN^, keZ (2.10) 

(in the confining phase) since under the large aperio- 
dic transformations given by Eq. (|2.6|l q{n/T)q{Q) 
z~'^q{n/T)q{Q). This generates an exact selection rule in 
quenched QCD. The center symmetry is explicitly bro- 
ken by the presence of dynamical quarks and the choice 
of an order parameter for confinement is not obvious l4l| . 
As a consequence the selection rule implied by eq. (I2.10|) 
is no longer fulfilled. Nevertheless such selection rule be- 
comes relevant to chiral quark models in the large Nc 
limit and departures from it are found to be suppressed 
within chiral quark models in the large limit at low 
temperatures, due to the spontaneous breaking of chi- 
ral symmetry which generates heavier constituent quarks 
from light current quarks."^ This issue will be analyzed 
along this paper. 

C. Chiral quark models at finite temperature 

Chiral quark models have been used in the past to 
provide some semiquantitative understanding of hadronic 
features in the low energy domain. At zero temperature 
chiral quark models are invariant under global SU(7Vc) 
transformations. There has always been the question of 
how the corresponding constituent quark fields transform 
under local color transformations or whether a physi- 
cal gauge invariant definition can be attached to those 
fields ^31 ■ If we assume that they transform in the 
same way as bare quarks, it seems unavoidable to cou- 
ple gluons to the model in the standard way to main- 
tain gauge invariance as done in previous works (see 
e.g. Refs. 0j H)- These gluon effects are treated 
within perturbation theory at T = 0. This approxima- 
tion induces some sub-leading corrections in the calcu- 
lation of color singlet states where the effects of con- 
finement can be almost completely ignored for the low 
lying states This perturbative gluon dressing also 

complies with the interpretation that the whole quark 
model is defined at a low renormalization scale, from 
which QCD perturbative evolution to high energies pro- 
cesses can be successfully applied |46j. When going to 



^ In this formula (g(ra/T)q(0)) denotes contributions to the quark 
propagator including only paths which wind n times around the 
thermal cylinder. The average is for the quenched theory. 

^ We emphasize that our use of the approximate rule is in contrast 
to the so-called canonical ensemble description of QCD where, 
upon projection, triality is assumed to be exact even in t he pres- 
ence of dynamical quarks. See e.g. the discussion in |41| . 
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finite temperature, chiral quark models predict alre ady 
at the one loop level a chiral phase transition |22| 
at realistic temperatures. However, even at low temper- 
atures single quark states are excited what is obviously 
not very realistic for it means that the hot environment 
is in fact a hot plasma of quarks. On the other hand 
since the constituent quark mass is about a factor of 2 
larger than the pion mass, pion loops dominate at low 
temperatures Jgjj (for a review see e.g. Ref. p^ ). as ex- 
pected from chiral perturbation theory [T^ |2C| . In the 
present work we will deal with two chiral quark models, 
the Nambu-Jona-Lasinio (NJL) model 0)E,'49| where 
quarks are characterized by a constant constituent mass 
in the propagator due to the spontaneous breaking of chi- 
ral symmetry and the recently proposed spectral quark 
model (SQM) H^, El El where the notion of ana- 
lytic confinement is explicitly verified. For completeness 
we review briefly the corresponding effective action be- 
low. One common and attractive feature of chiral quark 
models is that there is a one-to-one relation to the large 
Nc expansion and the saddle point approximation of a 
given path integral both at zero and at finite tempera- 
ture. 



1. The NJL model 

The NJL Lagrangian as will be used in this paper reads 
in Minkowski space'' 

-Cnjl g(j^ - Mo)q 

Nf-l 

+ f E ((9^-9)' + (9^a»75g)') (2.11) 

a=0 

where q = {u,d, s, . . .) represents a quark spinor with 
Nc colors and Nf flavors. The A's are the Gell- 
Mann flavor matrices of the U{Nf) group and Mq = 
diag(rnu, rud, rus , . . .) stands for the current quark mass 
matrix. In the limiting case of vanishing current quark 
masses the classical NJL-action is invariant under the 
global U{Nf)ii (g) U{Nf)L group of transformations. Us- 
ing the standard bosonization procedure (54] it is conve- 
nient to introduce auxiliary bosonic fields {S, P, V, A) so 
that after formally integrating out the quarks one gets 
the effective action^ 

rNjL[5,P] = -^iVcTrlog(^D) 

- d^xtrf{S^ + P') . (2.12) 

We use Tr for the full functional trace, tr/ for the trace 
in fiavor space, and tr^ for the trace in color space. Here, 



We use Bjorken-Drell convention throughout the paper. 
^ Obviously at finite temperature the quark fields satisfy antiperi- 
odic boundary conditions whereas the bosonized fields obey pe- 
riodic boundary conditions. 



the Dirac operator is given by 

iB = - Mo - (S + i-fsP) ■ (2.13) 

The divergencies in Eq. H2.12|l from the Dirac deter- 
minant can be regularized in a chiral gauge invariant 
manner by means of the Pauli-Villars method, although 
the issue of regularization is of little relevance at finite 
temperature [22| for T A. This model is known 
not to confine and to produce a constituent quark mass 
M ^ SOOMeV due to the spontaneous breaking of chiral 
symmetry at zero temperature. The Goldstone bosons 
can be parameterized by taking 

S + iP = VU^VU (2.14) 

with U a unitary matrix (see Eq. H2.19|l ) with = E, 
and one can use that E = Af -t- with the scalar field 
fluctuation. The partition function for this model can be 
written as 

Znjl = y"-DJ7DEe*^"JL[^'^l . (2.15) 

By minimizing Fnjl one gets S = M, which generates 
the spontaneous breaking of chiral symmetry, and one 
obtains the gap equation 

where the Pauli-Villars regularization has been used. The 
Pauli-Villars regulators fulfill cq = 1, Aq = and the 
conditions q = 0, CiA| = 0, in order to ren- 
der finite the logarithmic and quadratic divergencies, re- 
spectively. In practice it is common to take two cut- 
offs in the coincidence limit Ai A2 = A and hence 
E,c,/(A2) = /(0)-/(A2) + AY(A2). 



2. The SQM model 

In the SQM the effective action reads 

rsQM[C/] - -iNc j dcjp(c^)Trlog(iD), (2.17) 

where the Dirac operator is given by 

iT> = i^-ujU'"--Mo (2.18) 

and p{uj) is the spectral function of a generalized 
Lehmann representation of the quark propagator with uj 
the spectral mass defined on a suitable contour of the 
complex plane [13 [H El HI . The use of certain spec- 
tral conditions guarantees finiteness of the action. The 
matrix U = = e^'^'^l f ( / is the pion weak decay con- 
stant in the chiral limit) is the flavor matrix representing 
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the pseudoscalar octet of mesons in the non-hnear repre- 
sentation, 









K+ 



















A judicious choice of the spectral function based on 
vector meson dominance generates a quark propagator 
with no-poles (analytic confinement). More details of 
the SQM at zero and finite temperature relevant for the 
present paper are further developed at Appendix^ The 
partition function for the SQM can be written as 

^SQM = / 7?C/e'rs«"[^l . (2.20) 



consequences of a spatially constant Polyakov loop. In 
this case the functional integration Dfl becomes a sim- 
ple integration over the gauge group dH.. The issue of 
locality is reconsidered in Section IVl CI 

III. UNNATURALNESS OF CHIRAL QUARK 
MODELS AT FINITE TEMPERATURE 

In this section we analyze the problem of chiral quark 
models at finite temperature, its interpretation in terms 
of thermal Boltzmann factors as well as the correspond- 
ing conflicts with Chiral Perturbation Theory at finite 
temperature. 

A. The problem 



D. The Polyakov-Chiral Quark Model 



As we will show in Sect. Illlj there is a conflict between 
large gauge invariance at finite temperature, reviewed in 
the previous Sects. Ill Al and lllBI and the standard chiral 
quark models presented in Sect. Ill CI The chiral quark 
model coupled to the Polyakov loop that will be moti- 
vated in Sect. II VI and analyzed in the rest of this paper 
synthesizes the solution and corresponds to simply make 
the replacement 



Oa Oa— IAa 



(2.21) 



in the Dirac operators, eq. H2.13|l and eq. H2.18|l . and 
integrating further over the A4 gluon field in a gauge in- 
variant manner [ssf yielding a generic partition function 
of the form 



(2.22) 



where DU is the Haar measure of the chiral flavor group 
SU{Nf )ii X SU{Nf)L and DQ the Haar measure of the 
color group SU{Nc), Tq is the effective gluon action 
whereas Fq stands for the quark effective action. If the 
gluonic measure is left out A4 — and = 1 we re- 
cover the original form of the corresponding chiral quark 
model, where there exists a one-to-one mapping between 
the loop expansion and the large Nc expansion both at 
zero and finite temperature. Equivalently one can make 
a saddle point approximation and corrections thereof. In 
the presence of the Polyakov loop such a correspondence 
does not hold, and we will proceed by a quark loop ex- 
pansion, i.e. a saddle point approximation in the bosonic 
field U, keeping the integration on the Polyakov loop fl. 
The work of Ref. Q corresponds to make also a saddle 
point approximation in Jl. In Section we stick to the 
one loop approximation and keep the group integration. 
This is the minimal way to comply with center symme- 
try at low temperatures. Although in principle Q{x) is 
a local variable, in what follows we will investigate the 



As already mentioned, chiral quark models at finite 
temperature have a problem since, even at low temper- 
atures, excited states with any number of quarks are in- 
volved, whether they can form a color singlet or not. 
This is hardly a new observation, the surprising thing 
is that nothing has been done about it so far, attribut- 
ing the failure to common diseases of the model, such as 
the lack of confinement. To illustrate this point in some 
more detail we will use a constituent quark model like 
the NJL model, where the quark propagator has a con- 
stant mass. To be specific, let us consider as an example 
the calculation of the quark condensate for a single fla- 
vor in constituent quark models with mass M. At finite 
temperature in the Matsubara formulation we have the 
standard rule 



dko 
'2^ 



F{k,ko) iT 



n— — 00 



F{k,iujn) 



(3.1) 



with ujn the fcrmionic Matsubara frequencies, a;„ = 
27rr(n-|- 1/2). For the discussion in this and forthcoming 
sections it is convenient to elaborate this formula a bit 
further. Using Poisson's summation formula 

F{m) ^ J2 da;i^(a;)e*2'^^" (3.2) 

m — — 00 n— — 00 

one gets the rule 

/ ^Fik,ko)^^ f: i^ir I ^Fik,^k,)e-'^^/^ . 

(3.3) 

In terms of the Fourier transform, one obtains for a finite 
temperature fermionic propagator starting and ending at 
the same point, 

00 

F{x;x)^ J2 i-'^TFi^^xo + in/T;x,xo). (3.4) 



n— — 00 



Note that the zero temperature contribution corresponds 
to the term n = in the sum. From a path integral 
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point of view, the zero temperature term comes from 
contractile closed paths whereas thermal contributions 
come from closed paths which wind n times around the 
thermal cylinder. For fermions, each winding picks up a 
— 1 factor. 

For the (single flavor) condensate we get^ 



(qq)* = -iN, J2 (-l)"trDirac^(a 



iCQ—in/T 



co—in/T 



(27r)4 P _ M2 

n— — oo " ^ ' 

{qq)^2^^^^Y.tiyiK,inM/T). (3.5) 



In writing the previous formula, finite cut-off corrections, 
appearing in the chiral quark models such as the NJL 
model at finite temperature have been neglected. This 
is not a bad approximation provided the temperature is 
low enough T ^ A (typically one has Awl GeV so even 
for T M K, 300 MeV the approximation works). At 
low temperatures we may use the asymptotic form of the 
Bessel function Hal 



to get for the leading contribution, 



(3.6) 



fMT\ 



3/2 



-MIT 



(3.7) 



This means a rather flat dependence on temperature for 
T M. (Numerically, the correction is about 1% for T « 
100 MeV for M = 300 MeV and {qq) « -(240 Me V)^). 
The strong attractive interaction which causes chiral dy- 
namical symmetry breaking is reduced at finite tem- 
perature and the energy is decreased by a decreasing 
constituent quark mass M*, eventually leading to a chi- 
ral phase transition |^ , the critical temperature is 
T K, 200 MeV.'' The coincidence of this number with lat- 
tice simulations has been considered a big success of chi- 
ral quark models and has triggered a lot of activity in the 
past (see e.g. Ref. and references therein). We show 
below that this apparent success might be completely 
accidental, as it does not incorporate basic physical re- 
quirements. 



B. Interpretation 

An interpretation of the previous formula for the con- 
densate is in terms of statistical Boltzmann factors. Us- 
ing the definition of the quark propagator in coordinate 
space 



Six) 



(27r)4 ^ - M 



= (i^ + M)A{x) (3.. 



with 
A{x 



d'^k e-''^--" _ KiiV-M'^x^) 



,(3.9) 



(27r)4 fc2 - M2 Aw^i y/-M^x^ 

at low temperature we get 

S{x,i/T)^ e-^^'^ . (3.10) 

Thus, for {qq)* and up to prefactors, we have the expo- 
nential suppression for a single quark propagator at low 
temperature. Using Eq. (|3.5|) and Eq. 1)3. 6|l the quark 
condensate can be written in terms of Boltzmann factors 
with a mass formula A'/„ — nM corresponding to any 
number of quark states. 

One might object against the previous interpretation 
by arguing that these factors only reflect in the Eu- 
clidean coordinate space the pole of the propagator in 
Minkowski momentum space, and hence that they are a 
natural consequence of the lack of confinement. While 
the former statement is true, in the sense that singulari- 
ties in Minkowski momentum space can be seen at large 
Euclidean coordinate values, the conclusion drawn from 
there is incorrect. As shown in Ref. (see Appendix^ 
quark propagators with no poles but cuts can also pro- 
duce a Boltzmann factor without prefactors, as it should 
be.^ To the same level of approximation, i.e. one quark 
loop, in the SQM we get (see Appendix 1X1 for details) 



Ml 
m 



= tanh(Ms/4r) 

= 1 - 2e-*'^^-/2T ^ 2e-^-^-/^ + • • • (3.11) 



where the "Boltzmann" constituent mass can be iden- 
tified with half the scalar meson mass M = Ms 1^2.? 
These calculations illustrate our main point and can be 
extended to any observables which are color singlets in 
the zero temperature limit; quark model calculations at 
finite temperature in the one loop approximation gene- 
rate all possible quark states. 



0,e-*^/^ + 0,,e-2^^/^ + --- (3.12) 



In what follows we use an asterisk as upperscript for finite tem- 
perature quantities, i. e. O* = Ot- 

The minimization can be written as the equation 
{qq)*{M*)/M* = {qq)(M)/M, so one has to know the 
mass dependence of the condensate at zero temperature. 



* Actually, the previous counter example shows that the lack of 
confinement has more to do with the presence of the exponential 
prefactors which are related to the available phase space. 

^ This relation together with the large Nc quark-hadron duality 
relation Ms = My discussed in Ref. \5?H yields M = My jl ~ 
385MeV, a reasonable value. 
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While there is no doubt that the leading term Oq has a 
Boltzmann factor corresponding to a single quark state, 
the term with mass 2M could in principle be equally a qq 
diquark state or a qq meson state. The latter possibility 
should be discarded, however. At one loop a qq pair 
can only come from the quark line going upwards and 
then downwards in imaginary time propagation. Since 
such a path does not wind around the thermal cylinder 
it is already counted in the zero temperature term. The 
qq contribution, instead, corresponds to the single quark 
line looping twice around the thermal cylinder and is a 
proper thermal contribution. This is confirmed below. 
These Boltzmann factors control the whole physics and 
temperature effects are sizeable for T M . 



C. Conflicts with ChPT 

Our observation on the Boltzmann factor is rather 
puzzling because it seems hard to understand how is it 
possible to generate non singlet states by just increasing 
the temperature. The reason has to do with the fact that 
the condensate itself is not invariant under Z(7Vc) trans- 
formations at finite temperature. For the example of the 
condensate we trivially obtain 

oo 

(99)* = E (-l)"(9>o)9(0)) (3.13) 

— XQ—in/T 
n— — 00 

i.e., the condensate at finite temperature can be written 
as a coherent sum of nonlocal quark condensates at zero 
temperature. If we make a gauge transformation of the 
central type, we get 



E 

n— — oc 



-z)"(g(xo)g(0)) 



XQ—in/T 



(3.14) 



i.e., the condensate can be decomposed as a sum of irre- 
ducible representations of a given triality n. Thus, the 
state with Boltzmann factor e""*^/^ is indeed a multi- 
quark state. 

This avoids the paradox, and suggests that in order 
to make a (centrally extended) gauge invariant definition 
of the condensate we could simply discard from the sum 
those terms which do not have zero triality, i.e. we would 
get 



00 

E 

n— — 00 



{~ir^^{q{x,)qm 



(3.15) 

iN^n/T 



This would generate as a first thermal correction a term 
with a Boltzmann factor corresponding to mass N^M (a 
baryon) which is obviously very much suppressed. Since 
a quark loop generates a dependence proportional to Nc 
we would obtain a Nce~^'^'^ dependence. 

Another problem now comes from comparison with the 
expectations of chiral perturbation theory at finite tem- 
perature In the chiral limit, i.e., for m^r ^ 2-kT <C 



47r/7r the leading thermal corrections to the quark con- 
densate for Nf = 2, for instance, are given by 



ChPT 



1 



8/1" 384/1 



(3.16) 



Thus, the finite temperature correction is A^c-suppressed 
as compared to the zero temperature value, since /^ 
scales as N^. This feature remains for finite pion mass, 
and is generic to any thermal correction in ChPT; the 
dominant contribution comes from quantum pionic fluc- 
tuations and not from quark thermal excitations. Al- 
though the previous formula predicts a lowering of the 
quark condensate, it cannot describe the chiral phase 
transition since ChPT assumes from the start a non van- 
ishing chiral condensate. In this sense, the scaling be- 
havior of the critical temperature with /^ and therefore 
with Nc suggested from direct extrapolation of the for- 
mula can only be regarded as an assumption. 

At this point we should remind that the mechanism by 
which the chiral symmetry is restored at finite temper- 
ature in standard chiral quark models in the one quark 
loop approximation is quite different from the trend de- 
duced from ChPT based mainly on pion loops. While 
in the first case it is due to populating states of the 
Dirac levels with the Fermi-Dirac thermal factor and a 
sudden decrease of the constituent quark mass gap 2M, 
in ChPT the "phase transition" is merely due to large 
quark-antiquark excitations with the lightest pion quan- 
tum numbers with a fixed gap (otherwise ChPT method 
cannot be applied) . These two pictures of the chiral sym- 
metry restoration are not dual to each other; the Nc be- 
havior of the critical temperature is different since in chi- 
ral quark models one has ~ M ~ N^ while in ChPT 
the extrapolated value of the "critical temperature" is 
Tc ^ 2\/2/7r ^ \/Ac. Quantum fluctuations have been 
included in chiral quark models at flnite temperature 
(for a review see e.g. Ref. jl^) and they are known to be 
1/A^c suppressed. Actually, the sub-leading 1/Nc contri- 
bution reproduces the first term of ChPT, eq. 1)3. 16|) . thus 
largely dominating at low temperatures. Taking into ac- 
count that ChPT by itself and more refined approaches 
incorporating meson resonance effects 57^ 58j provide a 
similar values of the "critical temperature" quite close to 
the lattice predictions for dynamical fermions and ex- 
trapolated to the chiral limit, one may wonder what is the 
meaning of the mean field quark chiral phase transition 
predicted in the past j23, |22, yj] and which has become 
a justification for chiral quark models at finite tempera- 
tures. These problems are also common to models where 
quarks and mesons are regarded as independent degrees 
of freedom. 

We will see in the rest of the paper how a convenient Nc 
suppression of quark thermal corrections arises naturally 
when a color source, the Polyakov loop, is coupled to the 
chiral quark model and subsequent projection onto color 
neutral states is carried out. In this scenario one would 
have a large transition temperature Tc ~ NcM due to 
quarks, i.e. no symmetry restoration due to filling in the 
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states above the Dirac levels in the absence of dynami- 
cal gluons and in the quenched approximation (Polyakov 
cooling). Non perturbative gluonic corrections modify 
this picture; they do predict instead a critical tempera- 
ture roughly equal the deconfinement phase transition, 
Tc = Td- Finally, pion loops are protected from addi- 
tional suppressions, so that the final result will be fully 
compatible with the ChPT behavior at low temperature. 



IV. COUPLING THE POLYAKOV LOOP IN 
CHIRAL QUARK MODELS 

A. General considerations 

As we have said, one can formally maintain gauge in- 
variance at zero temperature by coupling gluons to the 
model. In the spirit of the model these degrees of free- 
dom should be treated within perturbation theory, since 
the constituent quarks carry some information on non- 
perturbative gluon effects (see e.g. Ref. 0; Q| for ex- 
plicit calculations in the low energy limit). At finite 
temperature the situation is radically different; a per- 
turbative treatment of the Aq component of gluon field 
manifestly breaks gauge invariance (namely, under large 
gauge transformations). The consequences of treating 
such a coupling non-perturbatively in the case of a con- 
stant Aq field are straightforward and enlightening (see 
below for a discussion on the x-dependent case). 

Actually, in a more general context, the Polyakov loop 
appears naturally in any finite temperature calculation 
in the presence of minimally coupled vector fields within 
a derivative expansion or a heat-kernel expansion ap- 
proach. In this case, as shown in |^ [s^, a local one 
loop quantity, such as the effective Lagrangian or an ob- 
servable, takes the form 

£(x) =^tr[/„(r!(x))0„(x)] , (4.1) 

n 

where tr acts on all internal degrees of freedom, n labels 
all possible local gauge invariant operators On{x) (i.e. 
containing covariant derivatives), possibly with breaking 
of Lorentz symmetry down to rotational symmetry, and 
/„(ri(x)) are temperature dependent functions of the 
Polyakov loop which replace the numerical coefficients 
present in the zero temperature case. In this general con- 
text Vl{x) would be the local Polyakov loop of all mini- 
mally coupled fields. In particular, a chemical potential 
would give a contribution e^^'^ . Here we can see the ne- 
cessity of the presence of in (|4.1f) : being ^ a constant, 
it gives no contribution in the covariant derivative and so 



in On (a;), therefore the chemical potential can only act 
through the presence of the Polyakov loop in the expres- 
sion. This consideration also illustrates the breaking of 
gauge invariance in a perturbative treatment of fi: e''^"^ 
depends periodically on the chemical potential, with pe- 
riod 27rir, this is a consequence of the coupling of /i to 
the integer quantized particle (or rather charge) num- 
ber. Such periodicity would be spoiled in a perturbative 
treatment. Note that such periodicity is equivalent to 
one-valuedness of the functions /„ in 14.1|l . 

B. Coupling the Polyakov Loop 

Coming back to chiral quark models with gluonic 
Polyakov loops, in fact, the analogy with the chemical 
potential has been invoked in a recent proposal of K. 
Fukushima , which suggests coupling chiral quark 
models to the Polyakov loop at finite temperature in this 
way. Our own approach is similar, except that, as in 
(|4.1|) . we consider a local Polyakov loop Vl{x) coupled to 
the quarks. This is what comes out of explicit one loop 
calculations within a derivative expansion approach at fi- 
nite temperature [H^ El| ■ In those calculations there 
is a loop momenta integration at each given x, and the 
Polyakov loop appears minimally coupled, i.e., through 
the modified fcrmionic Matsubara frequencies, 

w„ = 27rr(n+ l/2 + i>), (4.2) 

which are shifted by the logarithm of the Polyakov loop 

n = e*^^'^ , (4.3) 

i.e. v{x) = Ai{x) / {2nT) . In our considerations, this 
is the only place where explicit dependence on color de- 
grees of freedom appears, so it is useful to think of v as 
the corresponding eigenvalues. The effect of such a shift 
corresponds to change Eq. (|3.4|l into 

oo 

F{x]x)-^ (-17(f))"#(f,a;o + m/T;f,a;o).(4.4) 

n— — oo 

The interpretation of this formula can be visualized in 
Fig. n in a one quark loop with any number of exter- 
nal fields at finite temperature and with a non-trivial 
Polyakov line, the quarks pick up a phase (—1) due 
to Fermi-Dirac statistics, and a non Abelian Aharonov- 
Bohm^^ factor SI each time the quarks wind around the 
compactified Euclidean time. The total contribution to 
the diagram is obtained by summing over all windings 
and tracing over color degrees of freedom. 



As noted below, in a model with vector mesons, there would 
be a corresponding flavor Polyakov loop. Such a contribution is 
expected to be much suppressed due to the large physical mass 
of the mesons. 



After our work was sent for publication Refs. |60(, l6lj appeared, 
extending the results of Fukushima. 

This is an electric type of phase and not the standard mag- 
netic one. The name is nonetheless appropriate since this electric 
phase was discussed first in the original AB paper. 
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FIG. 1: Typical one quark loop diagram with a non trivial 
Wilson line. For n windings around the (7(1) compactified 
imaginary time the quarks get a topological factor SI" in ad- 
dition to the Fermi-Dirac statistical factor ( — 1)". Wavy lines 
are external fields. The total contribution to the diagram is 
obtained by summing over all windings and tracing over color 
degrees of freedom. 



C. Dynamical Polyakov loop 

The above prescription gives the contribution for a 
given gluon field configuration, of which we have only 
retained the Polyakov loop.^'^ The next step is to inte- 
grate the gluons according to the QCD dynamics. This 
implies an average over the local Polyakov loop with some 
normalized weight pi^l] x)dQ. Here dQ is the Haar mea- 
sure of SU(A^c) and p{Q; x) the (temperature dependent) 
probability distribution of fl(x) over the gauge group. 
The emergence of the Haar measure of the integral rep- 
resentation of the Yang-Mills partition function was ex- 
plicitly shown in Ref. 55]. Due to gauge invariance, p(ri) 
will be invariant under similarity transformations, and 
hence it is really a function of the eigenvalues of Q. In 
this section we will mainly remain within a quenched ap- 
proximation and so the weight follows from a pure Yang- 
Mills dynamics, in particular the weight will be x inde- 
pendent, as we do not consider external fields coupled 
to the gluons. In Yang-Mills dynamics (in four dimen- 
sions and three colors) it its known to take place a first 
order transition from a center symmetric phase to a bro- 
ken symmetry or deconfining phase. Note that this is a 
rather peculiar phase transition where the symmetry is 
restored below the critical temperature, just the opposite 
as the standard case. Since the transition is discontinu- 
ous in observables such as the expectation value of the 
Polyakov loop, the probability distribution p{fl) will also 
be discontinuous as a function of the temperature at the 



In addition, gluons appear also perturbatively through the co- 
variant derivative. This will produce perturbative gluon ex- 
change contributions as in the zero temperature case. We will 
not consider those in this work. 

In Sections lVl Cl and lVlll we will discuss some implications about 
local corrections in the Polyakov loop and unquenched results, 
respectively. 



critical temperature. In the confining phase p{Q) will be 
invariant under Z,{Nc), p{^) — p{zQ). In the deconfining 
phase, such symmetry is spontaneously broken and one 
expects the Polyakov loop to concentrate around one of 
the elements of the center, at random. The coupling of 
dynamical quarks favors the perturbative value Q = 1 
{A4 = 0) as follows from computations of the effective 
potential at high temperature 64, 65, 66]. So in that 
phase we expect to have il concentrated near ft — 1, 
which would be equivalent to no Polyakov loop in the 
calculation. 

Actually, one does not need the full distribution of 
on SU(iVc), but only the marginal distribution of eigen- 
values. Denoting the Polyakov loop average by ( ), we 
have for a quark observable 

Cix) - ^(tre/„(f7)) tvOnix) . (4.5) 

n 

Consistently with gauge invariance, the functions fn{^) 
are just ordinary functions fn{z) evaluated a,t z = Q (e.g. 
e^) hence, if e"^^ , j = 1, . . . , Nc are the eigenvalues of 

(±trj{n)) = / d^pm^Y^fie^fn 

= J[^P(^)f(^"^) (4-6) 

with 

p{4>) := / dnp{n)—Y,2n5{4>-cj>,). (4.7) 

JSU(N,) ^^c 

Equivalently, all that it is needed is the set of momenta 
of the distribution, (trc(fi")). 

D. Group averaging 

At sufficiently low temperature in the quenched theory 
we can go further on the analytical side, since the distri- 
bution of the Polyakov loop becomes just the Haar mea- 
sure in this regime. As it will be discussed in section lVl Bl 
this fact is justified with results based on strong coupling 
expansions and in one massive gluon loop approximation. 
Actually, from eq. H6.18|l we find that in observables such 
as the quark condensate, the effect of p{^) being differ- 
ent from unity is almost negligible for all temperatures 
below the transition, implying that a Haar measure dis- 
tribution is an extremely good approximation in the con- 
fined phase. We elaborate further on gluonic corrections 
m section rvTBl 

The corresponding density of eigenvalues of the 

SU(A^c) group is given by [erLlesI 

^2..(5:^.)nie'^^-e*i^n|^^ (4-8) 

i—l i<j i—1 
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so of H4.7|l is simply 

= 1 - ^^^^ cos(7V,0) . (4.9) 

Using this result one can easily deduce the following use- 
ful formulas for the average over the SU(iVc) Haar mea- 
sure 

(N,, n = (4.10) 

(tr,(-^2)")su(«.) = J -1 : n = ±N, (4.11) 
I 0, otherwise (4.12) 

When this is inserted in, e.g., Eq. H4.5|l . one finds that the 
effect is not only to remove the triality breaking terms, 
as in Eq. H3.15|l . but additionally, the surviving thermal 
contributions are Nc suppressed as compared to the naive 
expectation. This solves the second problem noted in 
Section Iml 

E. Polyakov cooling mechanism 

In an insightful work, Fukushima 6] has modeled the 
coupling of the Polyakov loop to chiral quarks, with em- 
phasis in the description of the deconfining and chiral 
phase transitions (or rather, crossovers). The fact that 
the critical temperatures for both transitions are nearly 
equal, according to lattice calculations ,40j . finds a natu- 
ral explanation in that model. This follows from what 
we will call the Polyakov cooling mechanism, namely, the 
observation that, upon introduction of coupling with the 
Polyakov loop, any quark observable at temperature T 
(below Td) roughly corresponds to the same observable 
in the theory without Polyakov loop but at a lower tem- 
perature, of the order of T/Nc, as already noted in [6^ . 
This is a direct consequence of triality conservation. As 
discussed for Eqs. H3.14|l and l|3.15|l at the end of Sec- 
tion lllll Boltzmann weights e~^^'^ are suppressed in fa- 
vor of e^^<^*^/^. An extreme example of cooling would 
come from considering a U(l) gauge theory in a con- 
fined phase in such a way that f2 is a completely random 
phase coupled to the quark. This would be equivalent 
to a uniform average of the variable v in Eq. (|4.2() in 
the interval [0,1]. Clearly, such an average completely 
removes the discretization of the Matsubara frequencies 
and gives back the continuum frequency of the zero tem- 
perature theory. The same extreme cooling would ob- 
tain in a \]{Nc) gauge theory. In the SU(A^c) case the 
average is not so effective since the phases corresponding 
to each of the Nc colors are not changed independently, 
owing to the restriction det = 1 . The cooling mech- 
anism will be substantially modified in the unquenched 
theory, since sea quark loops allow to create thermal (i.e., 
with n different from zero in e.g. Eq. ^AAl ) color sin- 
glet quark-antiquark pairs which propagate without any 
direct influence of the Polyakov loop. 

The way Polyakov cooling brings the chiral and de- 
confining critical points to coincide is as follows. In the 



chiral theory without Polyakov loop, the critical temper- 
ature of the chiral transition is such that T^^^ < Tjj yet 
NcT^^^ > To- Hence, in the theory with coupling to 
the Polyakov loop, one finds that for T < Td Polyakov 
cooling acts, (qq)* becomes roughly that of T/Nc which 
is below T^^^ and chiral symmetry is broken. On the 
other hand, for T > T^, Polyakov cooling no longer acts 
and quickly becomes unity, as in the theory without 
Polyakov loop at the same temperature; since T is above 
T^^^, chiral symmetry is restored. As a consequence the 
chiral transition is moved up and takes place at the same 
temperature as the deconfining transition, T^'^ w Td. 
This result is consistent with |70| where it is shown that, 
at least in the large Nc limit, confinement implies chiral 
symmetry breaking. 

We note a difference in our treatment of the Polyakov 
loop coupling and that in 0, namely, we use a local 
Polyakov loop subject to thermal and quantum fluctua- 
tions, as described by the distribution p(0; x)dn. This is 
in contrast with 6] where is global and does not fluctu- 
ate. Instead is determined through a mean field min- 
imization procedure plus a specific choice of the allowed 
values (orbit) of Q on the group manifold. In this way a 
model is obtained which is simple and at the same time 
can be used to address key issues of QCD at finite temper- 
ature. Nevertheless let us argue why such an approach 
needs to be improved. At sufficiently low temperature 
the model in Ref. Q for the gluon dynamics consist just 
of the invariant Haar measure on the gauge group, there- 
fore any group element is just as probable as any other. 
If one takes some coordinates on the group manifold and 
makes a maximization of the resulting probability den- 
sity, one is actually maximazing the Jacobian and the 
result will depend on the coordinates chosen. In the de- 
confined phase the local Polyakov loop is still subject to 
fluctuations (even in the thermodynamic limit). A differ- 
ent quantity is fi, the spatial average of the local loop.^^ 
This is a global object by construction. Both quantities, 
fl{x) and n, have the same expectation value, due to 
translational invariance, but does not fluctuate in the 
thermodynamic limit. The usual effective potential is so 
defined that its minimum gives the correct expectation 
value, and so fl, but it does not give information on the 
fiuct nations of ^l{x). 

In the confining phase of the quenched theory triality is 
preserved, hence, after gluon average, Eq. H4.4|l becomes 

F{x;x)^ (4.13) 

OO 

{{-^ix)r^')F{x, xo + inNc/T- f , x^) , 

n— — OQ 

which is the quenching invoked in Section ITTTI 



The quantity Q so defined does not fie on the group manifold, so 
some prescription should be devised to map it onto the group. 
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V. ONE QUARK LOOP RESULTS 



The calculations outlined above in Sect. lIVI can be rou- 
tinely applied to all observables. A more thorough and 
systematic study will be presented elsewhere. As an illus- 
tration we show here low temperature results (i.e. retain- 
ing only the Haar measure in the gluon averaging) for the 
quark condensate and the pion weak and electromagnetic 
anomalous decays for their relevance in chiral symmetry 
breaking, both for the NJL model as well as for the SQM 
at the one quark loop level. In Section FVI Al we discuss 
the structure of higher order corrections due to quark 
loops while in Section IVl BI dvnamical gluonic effects are 
considered. Corrections beyond the quenched approxi- 
mation will be explicitly computed in Section IV 111 In 
Ref. [tJ we compute the full chiral Lagrangian at finite 
temperature at the one quark loop level. 



A. Results for Constituent Quark Models 

To visualize the additional suppression we apply the 
previous result to the calculation of the condensate at 
finite temperature. At the one loop level we just make 
the substitution iVc(-l)" ^ trc((-l^)"). We get 



This yields 



iAM tre((-l])") 



n— — oo 



(27r)4 fc2 - M2 



Ki{NcM/T)- 



(5.2) 



The dots indicate higher gluonic or sea quark effects. Be- 
cause T is small we have further 



(qq)* - {qq)+4: 



/ MT 



3/2 



(5.3) 



When compared to the ChPT result Eq. H3.1(j|l we see 
that the Nc suppression of the constituent quark loop 
model is consistent with the expectations. 
For the pion weak decay constant we obtain 



(27r)4 [k^ - M2]2 



yielding 

f*2 Af2 

The TT — !■ 77 amplitude is given by 

p* = 



XQ—in/T 

(5.4) 



(5.5) 



X J2 ti-c((-f^)' 



^—ik-: 



(27r)4 [fc2 - M' 



213 



xo—in/T 

(5.6) 



Using the value obtained at zero temperature, ^^77 
l/47r2/7r, consistent with the anomaly, we get 



F* 9M 

^ = i-^Lk,{n,m/t) 



(5.7) 



This obviously complies again to the fact that the leading 
low temperature corrections should be encoded in pionic 
thermal excitations rather than quark excitations. 



XQ—in/T 

(5.1) 



B. Spectral Quark Model 

In the spectral quark model one averages with a given 
spectral function the previous result (|3.11() and including 
the Polyakov loop average we get (sec Appendix ^ for 
details) 



-N,Ms/2T 



im) 

For the pion weak decay constant we obtain 



(5.8) 



f*2 1 



T ' 



(5.9) 



£^2 
J TT 



and the ir'^ 77 amplitude is given by 



■ 7r77 



= 1 



77 



1 



12 



T 



NcMy 

T 



(5.10) 



VI. CORRECTIONS BEYOND ONE QUARK sponds to the quenched approximation within the model, 

LOOP 



In the previous sections we have restricted to the one 
quark loop approximation for observables. This corre- 
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a b c 

FIG. 2: Typical higher quark loop diagram for the quark condensate operator qq. Quark lines with independent momenta may 
wind n-times around the compactified Euclidean time, yielding a Fermi- Poly akov factor (—$!)". Triality conservation allows the 
internal quark-antiquark lines to wind with opposite signs only once, yielding an exponential suppression e~^*^/^ for diagram 
a). A similar suppression occurs for diagram b) if the quark-antiquark windings happen at any of the bubbles. Diagram c) 
corresponds to summing up all intermediate states with the same quantum numbers and can be interpreted as a meson line. 



and to some extent this provides an oversimplified pic- 
ture. In the present section we discuss the kind of cor- 
rections that we expect to this approximation. 

A. Higher Quark Loop Corrections 

Going beyond the one quark loop approximation may 
require tedious calculations (see e.g. Refs. for ex- 

plicit calculations in the standard NJL with no Polyakov 
loop). However, some general features based on Nc 
counting rules at finite temperature can be deduced as 
follows. Let us, for instance, consider the three loop di- 



agram of Fig. lOi) contribution to the quark condensate 
in the NJL model in terms of quark propagators. Writ- 
ing out for simplicity the Matsubara frequencies only we 
have 

Fig.|2t) = S{w^^^)(g)S{w^^'>) (6.1) 

where (g) means tensor product in the Dirac and inter- 
nal space sense. Using the Poisson's summation formula, 
Eq. H3.2|l . and going to Euchdean time space we get 



ni,n2,n3 
/•oo 



^5-jni+"2+n3^g-^A//T(|ni|-|-|n2|-l-|n3l) ("62) 



r 



For this diagram triality conservation implies, n\ -\- ni-\- 
na — kNc and the minimum argument of the exponent 
corresponds to take ni = n2 = — 0, which is the zero 
temperature contribution. The next thermal correction 
at low temperature is given by ni = , ^2 = ^ = 1 
so the 3 loop diagram of Fig. (j^i) is suppressed by a 
thermal factor e~^*^/^, to be compared to the one quark 
loop suppression g"^"^*^/^. A similar thermal suppres- 
sion is obtained by inserting the standard bubble sum- 
mation which can be coupled to meson quantum numbers 
transforming the argument of the exponent 2M Mqq. 
Obviously, this contribution becomes most important for 
the lightest pion state. Actually, the quark-meson dia- 



gram in Fig. ^Bp) looks as a two loop bosonized diagram 
as shown in Fig. (|3;). For such a bosonized diagram 
the previous argument becomes actually much simpler, 
since the number of loops equals the number quark prop- 
agators. The pion polarization operator, proportional to 
the pion propagator can then be taken at zero temper- 
ature, since the most important suppression comes from 
the quark lines not coupled to pion quantum numbers. 

For a bosonized diagram with L quark loops we have to 
consider L-fold Matsubara generalization of the previous 
one quark loop correction Eq. (|4.4|) . Actually, the anal- 
ysis becomes simpler in coordinate space. Regardless of 
the total number of quark propagators we may choose to 
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apply the Poisson's summation to L quark propagators. 
This can be seen by just using the formula 



and ns = • • • ri-L = 0. Assuming Nc = 3 we have 



16 



E 



n.m— — oo 



1/T 



dx4F{x4 + n/Tm/T) 



^ / dxiF{x4, + n/T) (6.3) 

1 — ^ ~oo 



and its multidimensional generalization both in the sum 
and in the integral sense. This effectively means that 
it is possible to remove as many Poisson summations as 
coordinate integrals appear in the expression. Using L = 
I - {V -I) and AV ^ E + 2I we also have 



n 



E \{^^^{-^rS{x,,U + in,/T). 

l,...,nL 

(6.4) 

Actually, this rule does not depend on the precise form of 
the quark interaction. At low temperatures, each quark 
line with an independent Poisson index generates a con- 
stituent quark mass suppression. Thus, the contribution 
to an observable can schematically be decomposed as fol- 
lows 



+ \nL\)M/T 



L ni,...,nL 



(6.5) 

Triality conservation of the measure ^l ^ zft at this level 
implies 



ni + ■ ■ ■ + riL ^ Nek 



(6.6) 



with k an integer. The dominant term in the previous 
expansion is the one for which rii = . . . = = with 
any arbitrary number of quark loops L and corresponds 
to the zero temperature contribution. One also sees that 
for L = 1 we only have contributions from m = kNc, 
which give the correction ^-^cM/t ^ hence reproducing 
the results of Sect. According to Eq. (|6.5|) . we can 
organize the thermal expansion at finite but low temper- 
atures. The most important contributions comes from 
minimizing X^iLi l^-d subjected to the triality constraint, 
Eq. H6.6|l . At finite T and for A^'c > 3 we have the leading 
temperature dependent contribution is given by i > 2 
and rii = — n2 ~ 1 with n-^ = ■ ■ ■ = ul = 0, which gives 
a factor e^^^^/^ and corresponds to a qq singlet meson 
state. This contribution has an additional l/Nc power 
suppression, as compared to the zero temperature contri- 
bution. For Nc — 3 the next term in the expansion would 
correspond to L > 3 and ni = n2 = ns = 1 and yields a 
finite temperature suppression g-'^cM/T ^ ]\[^ > 5 we 
would instead get L > A and ni = —n2 = = = 1 



Z- Lp-2M/T 

Nc 



N 



(6.7) 
(6.8) 
(6.9) 
(6.10) 

'}_^-{2Nm+NbN^)M/T -^-^n 



^-{2+Na)M/T 



\qq)«M(qqq)«B ^ ^ 

Obviously, for Nc = 3 the meson loop contribution dom- 
inates over the baryon loop contribution. The previous 
argument ignores completely the quark binding effects so 
we should actually consider the relevant meson mass m, 
thus in summary one would get 



(6.12) 



This is how quark-hadron duality works at finite tem- 
perature in chiral quark models. As we see contribu- 
tions of pion loops are the most important ones, even 
though they are 1/A^c suppressed. Higher meson states 
contribute next to the total observable at finite T. This 
is what one naively expects and it is rewarding to see 
that such a feature arises as a consequence of including 
the Polyakov loop into the chiral quark model and sub- 
sequent projecting onto the gauge invariant color singlet 
sector. 

Thus, at finite temperature there are the standard 
power like 1/Nce~'^'^/'^ suppression for meson loops ac- 
companied by an exponential suppression and a finite 
temperature exponential ^-^"^^l"^ for baryon loops. Ob- 
viously the most important contributions at large Nc or 
low T are those due to meson loops. We conclude from 
this discussion that thermal pion loops are protected. 

The previous discussion has concentrated on quark ob- 
servables. For an observable like the Polyakov loop one 
would have instead 



L ni,...,nL 

and 



+---+nL\„-{\ni\+--- + \nL\)M/T 



(6.13) 



1 + ni 



riL = Nek 



(6.14) 



The leading low temperature contribution (in this case 
there is no zero temperature term) is then of the type 
TT-i = —1, n2 = ■ ■ ■ = riL = Q, corresponding to a single 



In the case without Polyakov loop one would have Z^n^ (qq)iVm ~ 
— 1 — e-(2iVm+]V,)M/T ins^^eafj^ Qq t]-[g leading contributions are 
those corresponding to one quark state. 



14 



antiquark loop screening the charge of the test Polyakov 
loop. The leading term scales as and is controlled 

by the constituent quark mass. Unlike the quark con- 
densate case this behavior should remain unchanged by 
pionic loops. 



B. Gluonic Corrections 

Up to now we have chosen to represent the full dy- 
namical gluonic measure by a simple group integration. 
Unfortunately, we do not know at present any general 
argument supporting the idea that there is a low tem- 
perature exponential suppression of gluon degrees of free- 
dom, leaving only the Haar measure as the only remnant 
of gluon dynamics. However, results based on strong 
coupling expansions 1681 l72l and in one massive gluon 
loop approximation |1(tI . Il8l | do provide such a suppres- 
sion and indeed recent lattice findings confirm a striking 
universality in all group representations and favoring the 
simple group averaging dominance mechanism in gluody- 
namics below the phase transition ■ More specifically, 
one finds both from lattice calculations ^3 ^^'^ from the 
group measure that 



(t7c f7) = 



(6.15) 



in the confining phase for the Polyakov loop in the ad- 
joint representation. (In the group integration case, the 
previous formula follows from 1)7. 6|) below). We stress 
that this result is not a consequence of triality preserva- 
tion since O is invariant under 't Hooft transformations. 
The previous equation is equivalent to (|trcOp) = 1. We 
note in passing that in the mean field approximation 
(ItrcOp) vanishes instead, due to the absence of fluctua- 
tions. 

We analyze now the two above mentioned models. 



1. Strong coupling expansion 

The gluon potential at the leading order result of the 
strong coupling expansion, for Nc = 3, is taken as l6^l73| 

-iVain] = Vgincl^] -a^/T = -2{d- 1) c^'^'^/^jtrcf^f 

(6.16) 

with the string tension a = (425 MeV)^. At the mean 
field level V^iuc leads to a first order phase transition with 
the critical coupling 2{d~ l)e-'^"/^° = 0.5153. One can 
fix the deconfinement transition temperature as the em- 
irical value To = 270 Me V by choosing a'^ = 272 Me V 
The corresponding mass is mc — aa — 664 MeV. 
At low temperatures we may expand the exponential in 
powers of the gluon action. 



l + iTr 



1 



(6.17) 



which introduces an exponential suppression for e '"c/r^ 
For a treatment based on an average over the Polyakov 



loop, the normalized weight p{^)d^l suggested by the 
strong coupling expansion will be 

p(f^) = iVexp(2(d-l)e-'"<=/^|tref7|2) , (6.18) 

where TV is the normalization constant. Such distribution 
preserves exact triality. At low temperature p(r2) is close 
to unity and the distribution coincides with the Haar 
measure, hence fl is completely random with equal pro- 
bability to take any group value. At higher temperature 
p(r2) tends to favor concentration of near the central 
elements of the group, with equal probability. 

This provides the following mass formula for the Boltz- 
mann argument of the exponential (in the notation of 
subsection IVI All 



M = nNcMg + niMqq 



Imc 



(6.19) 



which clearly shows that the leading thermal contribution 
at low temperatures is, again, provided by pion thermal 
loops, corresponding to n — I — Q and to = 1 due to 
NcMq ^ TOG ^ Mgg = rriT^. Note that numerically, even 
the two pion contribution would be more important than 
gluonic corrections. 



2. One massive gluon loop approximation 



In a series of recent works [la, ll^ the equation of state 
has been deduced for a gas of massive gluons with a tem- 
perature dependent mass in the presence of the Polyakov 
loop, reproducing the lattice data quite accurately above 
the deconfinement phase transition. The vacuum energy 
density reads 



trc log 



1 



(6.20) 



where uik = \/k^ 



with Tog the gluons mass and f2 



and trc are the Polyakov loop and the color trace in the 
adjoint representation respectively. This expression was 
discussed with a temperature dependent mass in the de- 
confined phase given by the plugging the Debye screening 
mass tog(T) = Tg{T)\/2, which at the phase transition, 
T = Tc, takes the value tog(7;) = 1.2 - 1.3Tc. It is 
worth noticing that, if one assumes a constant value for 
the gluon mass below the phase transition one gets at 
low temperatures 



^1 9 r d^k 



,-nLOk/T 



where the identity 



tr,17" tr,17" - 1 



(6.21) 
(6.22) 



has been used. Using the asymptotic representation of 
the Bessel functions we see that, up to pref actors, a sim- 
ilar suppression of the sort described in the strong cou- 
pling limit. Sect. IVIB II takes place. 
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C. Local corrections in the Polyakov loop 

Up to now we have assumed a constant D, field in 
space in our calculations. Quite generally, however, the 
Polyakov loop depends both on the Euclidean time and 
the space coordinate, as it comes out of explicit one loop 
calculations within a derivative expansion approach at fi- 
nite temperature |6^ . In the Polyakov gauge the 
temporal dependence becomes simple, but there is still an 
unknown space coordinate dependence. In such a case, 
the previous rules have to be modified, since Polyakov 
loop insertions carry finite momentum, and the result 
depends on the ordering of these insertions. If we still 
assume that the Polyakov loop is the only color source in 
the problem, we are naturally lead to consider Polyakov 
loop correlation functions. In the confining phase we ex- 
pect a cluster decomposition property to hold for any pair 
of variables. A convenient model to account for Polyakov 
loop correlations is 



-cr\x~y\IT 



(6.23) 



with a the string tension. This includes the correct 
screening of the color charge at large distances due to 
confinement and is consistent with (|7.7|l for two Polyakov 
loops at the same point. Thus, very different values of the 
spatial coordinate are suppressed, and it makes sense to 
make a sort of local approximation within the correlation 
length, and expand correlation functions in gradients in 
that limited region of space. Effectively, this corresponds 
to replace the volume to given confinement domain, by 
means of the rule 



Y 
T 



(fa 



(6.24) 



In Ref. |7lj | we will see explicitly that when computing the 
low energy chiral Lagrangian by expanding the effective 
action in derivatives of the meson fields there appear also 
gradients of the Polyakov loop. Actually, since we cou- 
ple the coordinate dependent Polyakov loop effectively as 
a a;-dependent color chemical potential our approach re- 
sembles a non-abelian generalization of the local density 
approximation of many body physics in nuclear physics 
and condensed matter systems, very much in the spirit 
of a density functional theory. 



VII. RESULTS BEYOND THE QUENCHED 
APPROXIMATION AT LOW TEMPERATURES 

A. General remarks 

The full Polyakov-Chiral quark model is given in 
Sect. Ill Dl bv Eg. (|2.22f) . Therefore any expectation value 
is defined as 



(7.1) 



with Fg given in (|6.16() and Tq [U, Jl] the quark con- 
tribution to the full action, given by H2.12(l in the NJL 
model and (|2.17|l for the SQM case. In the latter model 
the full quark contribution coincides with the fermion de- 
terminant, while in the NJL model there is an additional 
term arising from the bosonization procedure. 



jrQ[u,n] 



Det(iD)f2 exp 



4G 



d'^xtvfiM - Mo)' 



(7.2) 

(Note that here we have included in D the color degrees 
of freedom.) 

In this section we gather all our results to provide an 
estimate of the Polyakov loop expectation value at low 
temperatures as well as the quark condensate. This is 
particularly interesting since in the quenched approxi- 
mation {tic^) = 0, due to triality conservation. The 
fermion determinant does not conserve triality, but we 
show below that at low temperatures the violation is ex- 
ponentially suppressed, so that it is still a good quantum 
number, and the Polyakov loop can be used as an order 
parameter for center symmetry in the same way as the 
chiral condensate provides a measure of chiral symmetry 
restoration away from the chiral limit. 

In order to go beyond the quenched approximation, we 
will evaluate the fermion determinant in the presence of a 
slowly varying Polyakov loop following the techniques de- 
veloped in our previous work 59] . According to our dis- 
cussion of Sect. lVI O of local corrections, such an approx- 
imation makes sense in a confining region where there are 
very strong correlations between Polyakov loops. In the 
presence of the Polyakov loop the quark contribution can 
be generally written as 



(7.3) 



where C is the chiral Lagrangian as a function of the 
Polyakov loop which will be computed at finite tempera- 
ture in Ref. 71] in chiral quark models for non-vanishing 
meson fields. For our purposes here only the vacuum 
contribution with vanishing meson fields will be needed. 



B. SQM model 

In our case it is simpler to consider first the SQM. We 
have 



e^i^«[^'"l =Det(zD)o = e^^*/^, 



(7.4) 



where V is the three dimensional volume and —B* the 
vacuum energy density at finite temperature in the pres- 
ence of the Polyakov loop. The result for B* is quite 
simple and is listed in ljA19|) in Appendix ^ At low 
temperatures we may expand to get, 



^VB'/T ^ ^VB/T 



l-^e-^^/^^tr.(l7 + f^-i) 



(7.5) 
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with M = My 12 the constituent quark mass in the SQM 
and — B is the vacuum energy density at zero temper- 
ature, B = M^N^NfJim-K'^ = (0.2GeV)'* for three fla- 
vors (see Appendix The calculation of observables 
requires the group integration formula (7j|, 



ikOjl 



(7.6) 



whence one gets for the constant Polyakov loop case 



(7.7) 



Note that the effect of ignoring the Polyakov loop (i.e., 
setting i7 = 1) promotes this result by two orders in Nc- 
In this model the average over pion fields is trivial since 
the vacuum energy density does not depend on U at the 
one quark loop level. Neglecting momentarily the gluonic 
corrections Tq, using the previous formulas and (|7.1|) we 
get the leading order result 



L : 



1 BV 



-Mv/2T 



(7.8) 



Note that at this order the contribution from the denomi- 
nator is trivial. As expected, triality is not preserved due 
to the presence of dynamical quarks, but the relevant 
scale is the constituent quark mass. In addition, note 
that since B is proportional to Nc there is an extra 1/Nc 
suppression. So the Polyakov loop can be effectively used 
as an order parameter. Actually, our calculation suggests 
that a low temperature calculation of the Polyakov loop 
in full QCD might provide a method of extracting a gauge 
invariant constituent quark mass. Proceeding in a similar 
way from the expression of the quark condensate (jAlSp 
we get the leading order contribution 



1 



2BV 



^-(Mv+Ms)/2T _^ , 



(7.9) 



It is noteworthy that the thermal correction scales as 
1/A^c {B scales as iVc), as in the ChPT case. This again is 
not just a consequence of triality, but requires the proper 
integration over the Polyakov loop manifold. The pres- 
ence of the (infinite) four- volume factor V/T has to do 
with our assumption on a constant Polyakov loop. As 
we have argued in Section FVI CI one has indeed a local 
Polyakov loop and the volume should be replaced accord- 
ing to the rule in Eq. H6.24|l by an effective confinement- 
domain volume. 

The first gluonic correction contributes in L as 
^-{Mv+2mG)/2T ^ quark condensate as 

^-(Mv+Ms+2mG)/2T 



C. NJL model 

The previous computation can also be considered 
within the NJL model. In this model the fermion de- 
terminant can be obtained by means of a derivative ex- 
pansion [s^le^. The result will be presented in Ref. 0|. 
Retaining only the vacuum contribution, which coincides 
with the result given in Eq. (3) of 0, we have 

Det(iD)o = exp [i J d^x {Cq{T = 0) + Cq{n, T)) 

where Cq{T — 0) is the zero temperature contribution. 
At low temperature, the thermal correction reads 



27r3 



Using the volume rule J d'^x Cq 



(7.11) 

iV/T)Cq, expand- 



ing Eq. (|7.10|l in powers of Cq{fl, T) and considering the 
group integration formula as above, we get the leading 
order result 



Nr. T 



27r3 



(7.12) 



(Since the NJL bosonization term in H7.2|) cancels in the 
calculation of observables, it needs not be included in 
this calculation. Also the gluonic corrections have been 
omitted. Their effect is discussed below). 

For the quark condensate we take into account the re- 
sult similar to Eq. H3.5|l but replacing (—1)" with (— il)", 
corresponding to the quark condensate for fixed Polyakov 
loop. Thus, including the leading fermion determinant 
contribution, using (|7.6|) . and taking into account that 
(trc^l) = (trcfi"^), we get for the single flavor conden- 
sate 



NfV 



(MT)3e-2AVT 



(7.13) 



Note that the Nf factor comes from the fermion deter- 
minant. As in the spectral quark model, the first gluonic 
correction contributes in L with q-(m+"^g)/t ^ ^^^^^ ^j^g 
quark condensate with g,-i^^^+"^G)/T ^ 

As we see, beyond the quenched approximation the 
Polyakov cooling persists although is a bit less effective 
as in the quenched case, and for instance the tempera- 
ture dependence of the low energy constants of the tree 

level chiral effective Lagrangian becomes L* — Li '^"^ 
e-Mv/T 



^"^ For the expectation value of a local observable 0{x), points out- 
side the volume V are not correlated and their contribution ap- 
proximately cancels in numerator and denominator. 



Actually we find a negative value for the SQM and positive for 
the NJL model. While based on color-charge conjugation sym- 
metry it can rigorously be shown that L must be real no proof 
exists to our knowledge that L > at any temperature, although 
lattice data |74| favor the positive case. 
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Finally, on top of this one must include higher quark 
loops, or equivalently mesonic excitations, from which 
the pions are the dominant ones. They yield exactly the 
results of ChPT '2^ for the chiral condensate {qq) and 
for the would-be Goldstone bosons, pions dominate at 
low temperatures. Thus, we see that when suitably cou- 
pled to chiral quark models the Polyakov loop provides 
a quite natural explanation of results found long ago on 
purely hadronic grounds ^ direct consequence of 

the genuinely non-perturbative finite temperature glu- 
onic effects. The expected leading correction effect on 
the Polyakov loop is also an additional exponential sup- 
pression ©(e"™"/^). 



VIII. IMPLICATIONS FOR THE PHASE 
TRANSITION 

The inclusion of the Polyakov loop has the consequence 
that one changes the one quark state Boltzmann fac- 
tor NcC^^^^'^ into (tYcfl) at low temperatures. In the 
quenched approximation one has (trj,f2) = 0, whereas the 
first non- vanishing contribution stemming from the Dirac 
sea behaves as (trcfi) ^ e~*^/^ due to the explicit break- 
ing of the center symmetry induced by the fermion deter- 
minant. Likewise, for the quark condensate (qq) the finite 
temperature correction changes iVcC"*^/"^ ^-im/t 
ter the Polyakov loop integration is considered. Tak- 
ing into account the large number of approximations and 
possible sources of corrections it is difficult to assess the 
accuracy of these Polyakov Chiral Quark Models, in s pite 
of the phenomenological success achieved in Refs. [ollSOjl 
within the mean field approach. Nevertheless, it is tempt- 
ing to see how these results may be modified not only 
at low temperatures but also in the region around the 
phase transition when the proper quantum and local na- 
ture of the Polyakov loop is considered. This requires go- 
ing beyond low temperature truncations like H7.12|l and 
(|7.13|) . Clearly a proper description would demand a 
good knowledge of the Polyakov loop distribution as a 
function of the temperature. Unfortunately, such a dis- 
tribution is poorly known and lattice simulations are not 
designed to extract it, since a subtle renormalization is- 
sue is involved 0, l^^. As a first step to investigate 
the phase transition in the Polyakov chiral quark model 
beyond the mean field approximation we just take the 
strong coupling model for the gluonic action of H6.18|l . 
Due to the rather large exponential suppression, this 
ansatz has the virtue of reducing to the Haar measure 
in the low temperature regime, and as a consequence the 
vanishing of the adjoint Polyakov loop expectation value 
observed in lattice calculations follows. In our view 
this is a compelling reason to go beyond mean field by in- 
tegrating over Polyakov loops. However, such a distribu- 
tion preserves center symmetry and would not generate 
a phase transition per se in gluodynamics. This is unlike 
the mean field approximation where the action is mini- 
mized by center symmetry breaking configurations. As 
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FIG. 3: Temperature dependence of the chiral condensate {qq) 
and Polyakov loop expectation value L = {tVcO,)/Nc in rela- 
tive units. The standard result for {qq)* corresponds to the 
pure NJL model uncoupled to the Polyakov loop. The result 
of L for gluodynamics within the strong coupling expansion 
is also displayed. We compare the mean field approach of 
Ref. where the Polyakov loop is classical and coupled to 
the quarks, with the integration over the Polyakov loop fi. 



discussed before a side product of this approximation is 
to miss the fluctuations and also to introduce an explicit 
coordinate dependence in the gauge group. In our model 
the breaking of the center symmetry is attributed only 
to quarks. As we will see this explicit breaking is rather 
large precisely due to the simultaneous restoration of the 
chiral symmetry, since the constituent quark mass drops 
to zero. The qualitative agreement with lattice calcula- 
tions in full QCD suggests that an important part of the 
physics has been retained by the model, leaving room for 
improvement in the Polyakov loop distribution. 

We will present calculations only for the NJL model. 
In practice, we use (|7.1|l . where the fermion determi- 
nant corresponds to Eq. (3) of |^ plus the volume rule 
(|6.24|) . The Polyakov loop integration is carried out nu- 
merically. Due to gauge invariance the Polyakov loop 
dependence is through its eigenvalues, and thus one may 
use the marginal distribution of eigenvalues H4.8|l , which 
for Nc = 3 amounts to two independent integration vari- 
ables. Full details are given in appendix IbI 

In Fig. O we show the effect on both the chiral con- 
densate {qq) and Polyakov loop expectation value L = 
{tTc^)/Nc within several schemes. In all cases we al- 
ways minimize with respect to the quark mass and use 
p(ri) in (|6.18() for all temperatures. We compare the 
standard NJL model with no Polyakov loop with the 
mean field calculation of Ref. 6] , which corresponds to 
minimize the vacuum energy as a function of the con- 
stituent mass and a given choice of the Polyakov loop 
matrix. We also compare with the result one obtains 
by integrating in the Polyakov loop instead and mini- 
mizing with respect to the quark mass afterwards. We 
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FIG. 4: Temperature dependence of d{qq)* /dT and dL/dT in 
the NJL model when the integration over the Polyakov loop Q 
is carried out. 

work in these calculations with the NJL model with 2- 
flavor, Nf — 2, and consider for the current quark mass 
matrix Mq = diag(m„,md) the isospin-symmetric limit 
with rriu — rud = ruq = 5.5 MeV. The zero temper- 
ature part of the effective action of Eq. (|2.12|) is reg- 
ulated by the Pauli-Villars method, with the cut-off 
Apy = 828 MeV, corresponding to a constituent quark 
mass M = 300 MeV. The coupling is G =13_13GeV-2, 
which is obtained from the gap equation (|2.16|l . These 
parameters reproduce the empirical values of the pion 
weak-decay constant and the quark condensate at zero 
temperature. Aspects of locality have been considered 
in the treatment of the NJL model with the integra- 
tion in the Polyakov loop, by introducing the volume 
rule (|6.24l) . where the string tension has been fixed to 
its zero temperature value a — (425 MeV) ^. It is also 
displayed in the figure the expectation value L in glu- 
odynamics within the model of Eq. Ht).16|l in the mean 
field approximation, which leads to a first order phase 
transition at To = 270 MeV. As we see the net effect of 
the Polyakov loop integration is to displace the transition 
temperature to somewhat higher values. So, the method 
based on the integration provides an effective cooling at 
higher temperatures for fixed parameters. As we can see 
m Fig. m the crossover transitions for the chiral conden- 
sate {qq} and for the Polyakov loop expectation value L 
coincide at the value Tc ~ 256 MeV. 

We have checked that a temperature dependence of the 
string tension may accommodate the unquenched lattice 
results [t^I, as we can see in Fig. |S| This provides a 
range of string tensions a — 0.181 ± 0.085 GeV^ with 
somehow account for an estimate of the uncertainty in 
the present model. In Fig. O the error band associated 
to such an uncertainty reflects a critical temperature of 
about Td = 250 ± 50 MeV. This is compatible with the 
large rescaling advocated in Ref. [b^. At present, and 
taking into account the many possible sources of correc- 



FIG. 5: Temperature dependence of the chiral condensate 
(qq) and Polyakov loop expectation value L = {trcO,)/Nc 
in relative units, in the NJL model when the integration 
over the Polyakov loop fl is carried out. The error bands 
are associated to an uncertainty in the string tension of 
a — 0.181 ± 0.085 GeV^. We compare with lattice data cor- 
respond to 2-flavor QCD, taken from |7^ . 



tions to our calculations we do not see how more accu- 
rate predictions could reliably be made in the context of 
Polyakov-Chiral Quark Models. Nevertheless the semi- 
quantitative success indicates that essential features for 
the center symmetry breaking phase transition are en- 
capsulated by these models, and further attempts along 
these lines should be striven. Nevertheless, it should 
be reminded that although the breaking of the center 
symmetry in this model is only attributed to the pres- 
ence of quarks, one also has a contribution from gluons. 
In this regard let us mention that ignoring the expo- 
nentially suppressed gluon action (|6.18|l in the averag- 
ing has almost no effect below the phase transition and 
shifts up the transition temperature by about 30 MeV, a 
value within our error estimate. Given the importance 
of quarks in the phase transition one may wonder if the 
temperature dependent volume enhances the breaking of 
the center symmetry. In fact, the volume at the tran- 
sition temperature is roughly equal to gluon volume 
in (|6.16|l . At low temperatures the exponential suppres- 
sion dominates in the Polyakov loop expectation value 
where the volume appears as a harmless pref actor, see 
e.g. (|7.12l) . The effect of replacing the temperature de- 
pendent volume by a constant one can be seen in Fig. |S1 
Again changes are within our expected uncertainties. 

As we have argued the expectation value of the 
Polyakov loop is rather small at temperatures well below 
the phase transition. The difference between the mean 
field and the direct integration can be best quantified at 
the level of the fluctuations. While at the mean field 
level the probability of finding a given Polyakov loop 
would be a delta function, one expects a spreading of 
such probability due to quantum effects. For TVc = 3 the 
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Polyakov loop contains two independent variables, which 
correspond to gluon fields in temperature units. 

n = diag (e*"*i , 6'"*=, e-''('^i+'^^)) (8.1) 

The joint distribution ^(01, can be factorized as a 
product of the purely gluonic and the quark determinant 
contributions (see appendix IB|| 

P{(f>l,<l>2) = PG (</>!, 02 )PQ (01, (/'2) (8.2) 

echoing the effective action displayed in Eq. H2.22|l in Eu- 
clidean space. Note that /9(0i,02) is not normalized to 
unity, instead its integral gives the full partition function 
(see appendix IB|| . As noted in Sect. IIV CI bv gauge in- 
variance the distribution is invariant under permutation 
of the three angles 0i, 02 and 03 = —01 — 02- The use of 
such a symmetry is that the trace of any arbitrary func- 
tion of the Polyakov loop f{il) (a one-body operator) can 
be averaged over the group by integrating out one angle, 
Eq. H4.6|) . Thus one obtains an equivalent one-body dis- 
tribution as 

p(0)(x^/ d0'pG(0,0')PQ(0>')- (8.3) 

It is interesting to compare how this distribution evolves 
across the phase transition, and to look for the effects 
generated explicitly by the fermion determinant. In 
Fig. we present such a comparison. Below the phase 
transition, and as already advanced in Sect. IIV CI the 
weighting function presents three maxima at equidistant 
values, as required by the center symmetry. In this case 
the quark determinant plays a negligible role, although a 
tiny, indeed exponentially small, center symmetry break- 
ing can be observed. As we see there appears an inter- 
esting concentration of angles in the region around the 
origin as the phase transition takes place. The quarks are 
very effective suppressing contributions not near 17 = 1. 
As a consequence the lack of the spontaneous breaking of 
the center symmetry in (I6.18|l becomes not very relevant 
for temperatures above the transition. 

A further trace of fluctuations can be seen by consider- 
ing higher group representations of the Polyakov loop. In 
Fig. [3we also show the expectation value of the Polyakov 
loop in the adjoint representation, (tic^) /{N^ — 1). Ac- 
cording to the lattice results of the matrix model in 
Ref. hj one has a vanishing expectation below the phase 
transition. As we have argued above, this feature is not 
preserved at the mean field level, where a non- vanishing 
value —1/{N^ — 1) is obtained instead (see Eq. (|6.22|) for 
the case n — 1). Considering the Polyakov loop inte- 
gration, as we do, complies with the lattice expectations 
and indicates that further developments should consider 
these constraints. The full fluctuation of the Polyakov 
loop, S, is defined by 

= (l + (t'^r,f])-(tr,r!)2)/iV|. (8.4) 



The fluctuation is also shown in Fig. |7| S goes to zero in 
the large T regime, and this is compatible with the fact 
that the one-body distribution yo(0) tends to concentrate 
near = as the temperature increases. In the second 
equality of Eq. (|8.4(l we have used the identity (|6.22l) with 
n = I. 



IX. CONCLUSIONS 

In the present work we have discussed how the problem 
of conventional chiral quark models at finite temperature 
may be overcome by introducing the Polyakov loop. In 
order to maintain gauge invariance at finite temperature 
some non-perturbative explicit gluonic degrees of free- 
dom must be kept. In practice, and in particular gauges 
such as the Polyakov gauge, the approach corresponds 
to treat the Aq component of the gluon field as a color 
dependent chemical potential in the quark propagator. 
This introduces, however, a color source which generates 
any possible color non-singlet states, calling for a projec- 
tion onto the physical color singlet states, or equivalently 
evaluating the path integral over the Aq field in a gauge 
invariant fashion. As such, the average includes both 
the gluon action and the quark determinant. Models for 
the gluonic part have been discussed on the light of pure 
gluodynamics results on the lattice. The net result is 
that, contrary to standard chiral quark model calcula- 
tions at finite temperature, no single quark excitations 
are allowed in physical observables. More generally, the 
leading thermal corrections at the one quark loop level 
start only at temperatures near the deconfinement tran- 
sition. Given the fact that this strong suppression effect 
is triggered by a group averaging of Polyakov loops we 
have named this effect Polyakov cooling of the quark ex- 
citations. Thus, and to a very good approximation, we 
do not expect any important finite temperature effect on 
quark observables below the deconfinement transition. In 
particular the chiral symmetry breaking transition can- 
not occur before the deconfinement transition. In such 
a situation the biggest change of observables such as the 
quark condensate should come from pseudoscalar loops 
at low temperatures and perhaps other higher meson 
resonances at intermediate temperatures. This is pre- 
cisely what one expects from ChPT or unitarized ap- 
proaches thereof which effectively include these loops on 
resonances. It is rewarding to see how in practice the 
apparent contradiction between chiral quark models and 
ChPT in the intermediate temperature region may be 
resolved by a judicious implementation of the Polyakov 
loop. 

The extrapolation of these ideas to the phase transi- 
tion is straightforward but more ingredients are needed. 
As an illustration we have investigated in a model the 
kind of effects one might expect from such a schematic 
Polyakov-Chiral Quark Model when both the quantum 
and local nature of the Polyakov loop are taken into ac- 
count. Several interesting features arise out of such an 
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FIG. 6: Temperature dependence of the one-angle Polyakov loop distribution p{(j)) of H8.8|l as a function of the angle. Dash- 
dotted lines: quenched result (pa is included, pQ is not). Center symmetry is preserved. Solid lines: unquenched result (both 
factors pG and pQ are included) in the NJL model. Center symmetry is explicitly broken. Three temperatures nearby the 
transition (255 MeV) are considered. For convenience all distributions have been normalized to unity. 
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FIG. 7: Temperature dependence of the Polyakov loop ex- 
pectation value in the fundamental, {tVcQ) /Nc, and adjoint, 
{tTc^)/{Nc — 1), representations and the total fluctuation S 
of the Polyakov loop, in the NJL model when the integration 
over the Polyakov loop is carried out. 



investigation. At low temperatures the Polyakov loop is 
suppressed exponentially in the constituent quark mass 
suggesting that eventually more accurate lattice measure- 
ments might provide a method to extract the constituent 
quark mass in a gauge invariant fashion. According to 
our analysis corrections to this leading behavior are pro- 
vided by pion loops. It would be extremely helpful to find 
a general theoretical setup where these chiral corrections 
might be reliably computed. Moreover, we find that the 
explicit breaking of the center symmetry due to dynami- 
cal quarks at low temperature is 1/iVc suppressed. This is 
a direct consequence of averaging over gauge field config- 
urations and confirms the current usage of the Polyakov 
loop as an order parameter in the unquenched case. On 
the light of the present findings one might conjecture that 



in the large Nc limit the Polyakov loop becomes a true 
order parameter of full QCD. 

Another feature we find is that the contribution of 
the gluon dynamics below the phase transition does not 
seem to be crucial. This is welcome since this is pre- 
cisely the region where least information can be deduced 
from lattice simulations besides the known preservation 
of the center symmetry. Nevertheless, it would rather 
interesting for our understanding of the low tempera- 
ture gluon dynamics to compute directly from the lattice 
the Polyakov loop probability distribution. From our re- 
sults we deduce that although the qualitative features 
observed in more simplified treatments are confirmed by 
calculations, one might expect large uncertainties in the 
determination of critical parameters, such as the critical 
temperature. Our estimate is To = 250 ± 50MeV for 
Nf = 2. Even given these large uncertainties, the very 
fact that a crossover between chiral symmetry restora- 
tion and center symmetry breaking takes place in the 
bulk part of the expected lattice QCD simulations with 
a minimal number of parameters is very encouraging and 
motivates that further studies along these lines should be 
pursued. 

Finally, a more intriguing aspect regards what kind 
of model independent information could be inferred out 
of these models, where quarks and Polyakov loops are 
coupled, in the regime around the phase transition. For 
instance, the low temperature behavior of the chiral con- 
densate can be described using Chiral Perturbation The- 
ory in terms of the zero temperature chiral condensate 
with no explicit reference to the underlying quark and 
gluonic degrees of freedom due to the dominance of pio- 
nic fiuctuations. Given the fact that the Polyakov loop is 
a gauge invariant object which vanishes at zero tempera- 
ture, it would be extremely helpful to isolate what physi- 
cal states could equivalently describe such an observable 
and what specific zero temperature QCD operators drive 
its low temperature behavior. 
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APPENDIX A: THE SPECTRAL QUARK 
MODEL AT FINITE TEMPERATURE 

In this Appendix we show how calculations for the 
spectral quark model introduced in Refs. 5(1 i51i i52l 
can be extended at finite temperature. 



1. The spectral quark model at zero temperature 



In the spectral quark model l5Ct |51| the quark propa- 
gator is written in the generalized Lehmann form 



S{k) 



c ft ~ u) 
Z{P) 



ft-M{k^) 



(Al) 



where /o(w) is the, generally complex, spectral function 
and C denotes a contour in the complex lu plane. Mlk"^) 
is the self energy and Z{k'^) is the wave function renor- 
malization. As discussed already in Ref. ^51] the proper 
normalization and the conditions of finiteness of hadronic 
observables are achieved by requesting an infinite set of 
spectral conditions for the moments of the quark spectral 
function, p{uj), namely 

Po = J duipiuj) = 1, (A2) 

Pn = J dLJUj"p{uj) = 0, for n = 1,2,3,... (A3) 

Physical observables are proportional the zeroth and the 
inverse moments, 



= dLOUj-'^piuj), for n = 0,l,2,... (A4) 



P- 

as well as to the "log moments", 

p'^ = [ dLu\og{LU^)uj"p{Lj), forn = 1,2,3,... (A5) 



Obviously, when an observable is proportional to the 
dimensionless po = 1 the result does not depend ex- 
plicitly on the regularization. The spectral conditions 
(|A3p remove the dependence on a scale p in IjASp . thus 
guarantying the absence of any dimensional transmuta- 
tion. No standard requirement of positivity for the spec- 
tral strength, p{uj), is made. The spectral regularization 



is a physical regularization in the sense that it provides 
a high energy suppression in one quark loop amplitudes 
and will not be removed at the end of the calculation. 
Using the methods of Ref. 'bS] one finds at zero tempera- 
ture, T = 0, and trivial Polyakov loop ~ I, and in the 
chiral limit, = 0, the following results for the pion 
weak decay constant, /^r, the single flavor condensate 
(qq), the energy- momentum tensor 9'^'^ and the anoma- 
lous vr" 27 amplitude 



Nf{qq) 



i2n) 



iiNcNf 
iiNcNf 



d^k 
(2^ 
d^k 



B{k') 
[kf^k^Aik^) 



d*k 



NcUJ (2^) 



[k'A{k')]' 



(A6) 
(A7) 



(A9) 



%) = 



The vacuum energy density is defined by ey — 
—B with B the bag constant. In the meson dominance 
(MD) version of the SQM one obtains for the even an 
odd components of the spectral function 



pvii^) = 



1 1 1 



27rit^ (l-4w2/M2)5/2 

pi-u;)] 
Up', 



lipi^ 
1 



(AlO) 



(All) 



27riM4(l-4cj2/Af|)'^/2 

{My and Ms are the vector and scalar meson masses, 
respectively) and hence 

1 



Thus, 



Aik') 
Bik') 



1 



(l-4fc2/M2)5/2 

487r2(g(j) 



M|iVc(l-4fc2/M|)5/2' 



247r2 



tv = -B = -'- 



N 



(A12) 

(A13) 
(A14) 



1927r2 

For three fiavors one has B ~ (0.2GeV)'* for My 
770MeV. 



2. The spectral quark model at finite temperature 
and arbitrary Polyakov loop 

We introduce finite temperature and arbitrary 
Polyakov loop by using the rule 

/rlh °° 
-^F{k,ko)^iT J2 nl^^n) (A15) 
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with LOn the fermionic Matsubara frequencies, tD„ = 
2TTT{n+l/2+i)), shifted by the logarithm of the Polyakov 
loop il — e'^^^'^ . In the meson dominance model, we use 
the momentum space representation, evaluate first the 
three dimensional integrals and finally the sums over the 
Matsubara frequencies. In practice all sums appearing 
are of the form 



C30 



i-i 



[i-iy. V dA'p^ 

The basic sum is given by 

1 sinh(Af/T) 



SiiM,T) 



2M cos(27rz>) + cosh(M/r) 



Si{M, r).(A16) 



(A17) 



Using the relations in [5l| and the previous formulas, we 
get 



Nr. 



sinh(M5/2T) 



cos(27ri>) + cosh(Ms/2T) 



(A18) 



and, for the vacuum energy density 



B* 



-tr. 



sinh(My/2T) 



cos(27r;>) + cosh(Afy/2r) 



(A19) 



Note that the relative temperature dependence of the 
bag constant and the quark condensate are alike in the 
present model, if we also interchange the vector and 
scalar masses. The T and O dependence of is 



f*2 



Tsinh(My/r) - Mv - cos(27rj>) [My cosh(My/2T) - 2rsinh(Afy /2r)] 
Nr " 2T{cos{2TTi>)+cosh{Mv/2T)f 

I 



(A20) 



To compute the group averages we observe that for a 
sum of the form of Eq. I|A17|I we can make z = e^'^ and 
t = e^*^/-^. Expanding in powers of t we get 



1/t - t 



z + l/z + t+ l/t 

OO 



(A21) 



and applying Eq. (|4.1U|) - (|4.12|) for the average over the 
SU(A'c) Haar measure one has 



A< 



1 - —t^^ 

Nc 



(A22) 



Thus, undoing the change of variables we get 



1 



-trc 



sinh(M/r) \ 



1 



-NaM/T 



Nc ''cos(27rz)) +cosh(A//T)/ N^ 

(A23) 

Beyond the quenched approximation, integrals may be 
computed at low temperatures using Eq. H7.7|l and gen- 
eralizations thereof. 



APPENDIX B: DETAILS ON NUMERICAL 
GROUP INTEGRATION 

In this Appendix we give details relative to the calcula- 
tion presented in Sec. I Villi for the Nambu-Jona-Lasinio 
model. 



The chiral condensate is obtained from maximization 
of the partition function Z, Eq. (|2.22l) . with respect to 
the constituent quark mass. To compute Z we need to 
carry out first the color group integration. To this end we 
consider the Polyakov gauge and parameterize the SU(3) 
Polyakov loop matrix as in Eq. (|8.1|) . The expression is 



2n 2tt 



Pg('/>i,02)pq('/>i,02) , 



where we have separated a gluonic distribution 



2n 2tt 



and a fermionic one 



(Bl) 



(B2) 



(B3) 



Note that in Tq [Q] we do not consider any dependence in 
the mesonic U fields, because we only retain the vacuum 
contribution, pa contains the Haar measure associated 
with the SU(3) group integration, as well as the gluonic 
corrections given in Eq. H6.16|l . i.e. 



Pg(01>2 



and 



1 
6 

X exp 



27 - Itrcf^l"* + 8Re((trcr2)^) - 18\trc^\^) 
(2(d-l)e-^'^/^|trcr!|2) , (B4) 



-i(l/'l+02) 



(B5) 
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FIG. 8: Comparison of observables computed using a con- 
stant correlation volume a^, — 272 MeV (dashed line) 
versus the temperature dependent volume V of Eq. 16.2411 
(solid line). 

The quark distribution pq(0i, 02 ) follows from the vac- 
uum contribution in Eqs. H7.2|l and H7.10|l . To obtain 
rQ[f7] we use, passing over to Euclidean space, 

tTgin] = -^(c,{T = Q) + C,{T,n) 

+^tr/(M-Afo)'). (B6) 
where the correlation volume V is given in 

(jn21I)- More- 
over, 

/<Pk 
(2^^' ' ^^^^ 



£,(r,0) = -2TNf J ^(^tr,log[l + f7e-^^/n 
+tic log [1 + e-^"/'^] ^ . (B8) 

Here Ek = VF+M^ is the energy of quasi-quarks and 
the constituent quark mass is M = niq — G{qq). For the 
zero temperature term we use the Pauli-Villars regular- 
ization. After momentum integration we get 

C,{T = Q) = -^^c.(A?+Af2)^log(A? + M2) 

(B9) 

For the temperature dependent part, after momentum 
integration and expanding the logarithm function, we ob- 
tain 

c,{T,n) = 4(AfT)2^L^i^2(nM/r) 

n=l 

X (tre(17") -f trc(f^"")) (BIO) 
In numerical calculations we check for convergence in 
theses sums. The expectation value of any observable is 
computed as in Eq. (|7.1|) . For any general function f{fl), 
the group averaging can be advantageously evaluated us- 
ing Eq. (|4.6|l . where the one-body distribution is given 
by Eq. H8.3|l . It proves convenient to evaluate the dou- 
ble integral as an iterated one, since with the exception 
of the adjoint Polyakov loop, observables depend on one 
angle only. A comparison of observables calculated us- 
ing a constant correlation volume a^, = 272 MeV @ 
and the temperature dependent volume V of Eq. (|6.24|l 
is shown in Fig. |S1 



[1] A. Casher, Phys. Lett. B83, 395 (1979). 

[2] F. Karsch, (1998), Prepared for 6th Hellenic School 

and Workshop on Elementary Particle Physics;, Corfu, 

Greece, 6-26 Sep 1998. 
[3] A. Gocksch and F. Neri, Phys. Rev. Lett. 50, 1099 

(1983). 

[4] F. Sannino, Phys. Rev. D66, 034013 (2002), 

hep-ph/0204174 . 
[5] A. Mocsy, F. Sannino and K. Tuominen, Phys. Rev. Lett. 

92, 182302 (2004), hep-ph/0308135 . 
[6] K. Fukushima, Phys. Lett. B591, 277 (2004), 

hep-ph/0310121 . 
[7] K. Fukushima, Phys. Rev. D68, 045004 (2003), 

hep-ph/0303225 . 
[8] A. Gocksch and M. Ogilvie, Phys. Rev. D31, 877 (1985). 
[9] P. N. Meisinger, T. R. Miller and M. C. Ogilvie, Nucl. 

Phys. Proc. Suppl. 129, 563 (2004), hep-lat/0309097 . 
[10] L. McLerran, Pramana 60, 575 (2003), 

hep-ph/0202025 . 
[11] U. W. Heinz, Nucl. Phys. A721, 30 (2003), 

,nucl-th/02 12004, . 



[12] C. P. Korthals Altes, ■hep-ph/0308229 

[13] H. Meyer, Nucl. Phys. B230, 31 (1984). 

[14] B. Svetitsky, Phys. Rept. 132, 1 (1986). 

[15] A. Dumitru and R. D. Pisarski, Phys. Lett. B525, 95 

(2002), hep-ph/0106176 . 
[16] P. N. Meisinger, T. R. Miller and M. C. Ogilvie, Phys. 

Rev. D65, 034009 (2002), hep-ph/0108009 . 
[17] A. Dumitru, Y. Hatta, J. Lenaghan, K. Orginos and 

R. D. Pisarski, Phys. Rev. D70, 034511 (2004), 
hep-th/0311223 . 
[18] P. N. Meisinger, M. C. Ogilvie and T. R. Miller, Phys. 

Lett. B585, 149 (2004), hep-ph/0312272 . 
[19] J. Gasser and H. Leutwyler, Phys. Lett. B184, 83 (1987). 
[20] P. Gerber and H. Leutwyler, Nucl. Phys. B321, 387 

(1989). 

[21] V. Bernard, U. G. Meissner and I. Zahed, Phys. Rev. 

D36, 819 (1987). 
[22] C. V. Christov, E. Ruiz Arriola and K. Goeke, Acta 

Phys. Polon. B22, 187 (1991). 
[23] W. Florkowski and W. Broniowski, Phys. Lett. B386, 

62 (1996), hep-ph/9605315. . 



24 



[49] 



M. Oertel, M. Buballa and J. Wambach, Phys. Atom. 
Nucl. 64, 698 (2001), hep-ph/0008131 . 
P. N. Meisinger and M. C. Ogilvie, Phys. Lett. B379, 
163 (1996), hep-lat/9512011 . 

O. Kaczmarek, F. Karsch, P. Petreczky and R Zantow, 

Phys. Lett. B543, 41 (2002), hep-lat/0207002 . 

E. Megias, E. Ruiz Arriola and L. L. Salcedo, JHEP 

0601, 073 (2006), hep-ph/0505215 . 

E. Megias, E. R. Arriola and L. L. Salcedo, 

hep-ph/0410053 , Mini- Workshop on Quark Dynamics: 

Bled 2004, Bled, Slovenia, 12-19 Jul 2004. 

E. Megias, E. Ruiz Arriola and L. L. Salcedo, 
[hep-ph/0411293 , 6th Conference on Quark Confinement 
and the Hadron Spectrum, Villasimius, Sardinia, Italy, 
21-25 Sep 2004. 

N. P. Landsman and C. G. van Weert, Phys. Rept. 145, 
141 (1987). 

C. Garcfa-Recio and L. L. Salcedo, Phys. Rev. D63, 
045016 (2001), hep-th/0007183 . 

A. N. Redlich, Phys. Rev. Lett. 52, 18 (1984). 

S. Deser, R. Jackiw and S. Templeton, Ann. Phys. 140, 

372 (1982). 

R. D. Pisarski, Phys. Rev. D35, 664 (1987). 

K. S. Babu, A. K. Das and P. Panigrahi, Phys. Rev. 

D36, 3725 (1987). 

G. Dunne, K.-Y. Lee and C.-h. Lu, Phys. Rev. Lett. 78, 
3434 (1997), hep-th/9612194 . 
L. L. Salcedo, Nucl. Phys. B549, 98 

hep-th/9802071 
L. L. Salcedo, 

hep-th/0201057 . 

G. 't Hooft, Nucl. Phys. B153, 141 (1979). 

F. Karsch, Lect. Notes Phys. 583, 209 
hep-lat/0106019 . 

K. Fukushima, Ann. Phys. 304, 72 

hep-ph/0204302 . 
M. Lavelle and D. McMuUan, Phys. Rept. 279, 1 (1997), 

hep-ph/9509344 . 

D. Espriu, E. de Rafael and J. Taron, Nucl. Phys. B345, 
22 (1990). 

J. Bijnens, C. Bruno and E. de Rafael, Nucl. Phys. B390, 
501 (1993), hep-ph/9206236 . 

L. Y. Glozman and D. O. Riska, Phys. Rept. 268, 263 
(1996), hep-ph/9505422 . 

E. Ruiz Arriola, Acta Phys. Polon. B33, 4443 (2002), 
hep-ph/0210007 . 

S. P. Klevansky, Rev. Mod. Phys. 64, 649 (1992). 

C. V. Christov et ai, Prog. Part. Nucl. Phys. 37, 91 

(1996), hep-ph/9604441 . 

R. Alkofer, H. Reinhardt and H. Weigel, Phys. Rept. 
265, 139 (1996), |hep-ph/9501213 . 



Phys. Lett. B530, 244 



(1999), 
(2002), 

(2002) , 

(2003) , 



[50 
[51 

[52; 



[53; 

[54' 
[55; 

[56; 

[57 

[58; 
[59; 

[60 
[61 
[62; 
[63 

[64 
[65^ 
[66 

[67 

[68 
[69 

[7o; 

[71 
[72; 
[73; 

[74; 



E. Ruiz Arriola, hep-ph/0107087 

E. Ruiz Arriola and W. Broniowski, Phys. Rev. D67, 
074021 (2003), hep-ph/0301202 . 

E. Ruiz Arriola and W. Broniowski, 'hep-ph/0310044' 
Light-Cone Workshop: Hadrons and Beyond (LC 03), 
5-9 Aug 2003, Durham, England. 

E. Megias, E. Ruiz Arriola, L. L. Salcedo and 
W. Broniowski, Phys. Rev. D70, 034031 (2004), 
hep-ph/0403139 . 

T. Eguchi, Phys. Rev. D14, 2755 (1976). 
H. Reinhardt, Mod. Phys. Lett. All, 2451 (1996), 
hep-th/9602047 . 

M. Abramowitz and L A. Stegun, Handbook of mathe- 
matical functions (Dover, New York, 1970). 
J. R. Pelaez, Phys. Rev. D59, 014002 (1999), 
hep-ph/9806539 . 

J. R. Pelaez, Phys. Rev. D66, 096007 (2002), 
hep-ph/0202265 . 

E. Megias, E. Ruiz Arriola and L. L. Salcedo, Phys. Lett. 
B563, 173 (2003), hep-th/0212237 . 

C. Ratti, M. A. Thaler and W. Weise, Phys. Rev. D73, 
014019 (2006), hep-ph/0506234 . 

S. K. Ghosh, T. K. Mukherjee, M. G. Mustafa and 
R. Ray, hep-ph/0603050 

E. Megias, E. Ruiz Arriola and L. L. Salcedo, Phys. Rev. 

D69, 116003 (2004), hep-ph/0312133 . 

E. Megias, hep-ph/0407052 9th Hadron Physics and 7th 

Relativistic Aspects of Nuclear Physics, Angra dos Reis, 

Brazil, 28 Mar-3 Apr 2004. 

N. Weiss, Phys. Rev. D24, 475 (1981). 

N. Weiss, Phys. Rev. D25, 2667 (1982). 

D. J. Gross, R. D. Pisarski and L. G. Yaffe, Rev. Mod. 
Phys. 53, 43 (1981). 

W. Miller, Symmetry groups and their applications (Aca- 
demic Press, New York, 1972). 
M. Gross, Phys. Lett. B132, 125 (1983). 
M. Oleszczuk and J. Polonyi, Annals Phys. 227, 76 
(1993). 

S. R. Coleman and E. Witten, Phys. Rev. Lett. 45, 100 
(1980). 

E. Megias, E. Ruiz Arriola and L. L. Salcedo, (2006), in 
preparation. 

J. Polonyi and K. Szlachanyi, Phys. Lett. BllO, 395 
(1982). 

M. Creutz, Quarks, Gluons and Lattices, Cambridge 
Monographs On Mathematical Physics, (Cambridge Uni- 
versity Press, Cambridge, UK, 1983). 
O. Kaczmarek and F. Zantow, Phys. Rev. D71, 114510 
(2005), hep-lat/05030171. 



